1959 


uc- 


1ESon 


yas ; 
ible | 
i- 


3 of 


ent 
ige- 
) 


the 
op- 

re- 
3sed. 


PER- 
As- 
r, 


fed 


Fermi 
8 are 


yhich | 


rons 
rmi 











PHYSICAL REVIEW 
LETTERS 








VoLuME 3 


JULY 15, 1959 


NUMBER 2 











CALCULATION OF PARTITION FUNCTIONS* 


J. HubbardT 
Department of Physics, University of California, Berkeley, California 
(Received June 15, 1959) 


It is the purpose of this Letter to describe 
briefly the manner in which a transformation 
due to Stratonovich' can be used to calculate 
approximately the grand partition functions of 
some many-body systems. 

Consider the grand partition function 


Z =Tr[exp(-8K -8V -vN)], 
where 


* * 
K= Lu; EM; "1, , V=$ Pe ais RI" my Ny» 


* 
N=) 1, n> B=1/KT, v=u/T, 


(T= temperature, » =chemical potential) and the 


7 and 7; are the creation and destruction opera- 
tors for the particles of the system; the n; * and 
hj may have commutation relations appropriate 
to either Bose-Einstein or Fermi-Dirac statis- 
tics, 

It should be noted that v; i: RL = Rl; ij necessar - 
ily, so that V can be regarded asa ’ quadratic 
form in the variables bij =" 1j- Bringing this 
quadratic form into normal form one has 
V=)),AaPa’, where the Aq are real and the py 
are linear combinations of the b5j- The grand 
partition function may now be rewritten 


Z=Tr[exp(-8K-vN-62),2,P,"))- (1) 


One can now transform this expression by 


making use of the identity 


f ” exp|- ax? - 2n”ax}dx =exp{a*}. (2) 


00 


This identity cannot be applied directly to (1) 
because of the noncommutativity of the operators 
involved. To overcome this difficulty one may 
make use of the ordering label technique of 
Feynman’ and introduce an ordering label s to 
write 


Z =Tr{exp[- 2 (K+HN )-L DAP 


a's @ as 


*I, (3) 


where )), has been written for f,°ds; the quanti- 
ties in (3) can now be handled as though they were 
c numbers, the ordering determined by s being 
restored when convenient. 

Applying the identity (2) to (3) one now easily 
finds 


z= f exp{- L[xDIax,., (4) 


where 


Lx 5 J=a 2 Xas + Aflx . 1 (5) 


and 


exp{- Affe I= T+ exp|-D 0K, + uN,) 


va 
29 F(-24) as? asl (6) 
The essential feature of the transformation 
(4)-(6) is that the operator appearing in the ex- 
ponent in (6) is only quadratic in the operators 
n;* » 1;, 80 that the evaluation of the indicated 
trace is essentially similar to the evaluation of 
the partition function of a system of noninteract- 
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ing particles moving in a field determined by the 
x... Equations (4) and (5) then express Z as a 
as 
weighted average of the quantity (6) over the 
fields [x5]. 
In some cases the integral (4) can be conven- 
iently evaluated by the method of steepest des- 


cents. One then finds 
InZ =-1).(x« °?- gflx °]-3 Infdet(Z,),], (7) 
aa a 


where the x,° are determined by the equation 


(8L/e% =0 (8) 


(it can be shown that the solution x,,,° of this 
equation is independent of s, hence the use of the 
abbreviated notation x,,° for xy,°), while L, is 
the matrix 


1 Ll 
Cos sys’ Qn oxox,” (9) 
as ys 
(),» means that the quantity is evaluated at x, 
2% ®, 

Evaluation of the formulas (7)-(9) shows that 
(8) represents an obvious generalization to finite 
temperature of the ordinary Hartree self-consist- 
ent-field equations, while the first two terms of 
(7) represent the energy in the Hartree approxi- 
mation and the last term of (7) is (at T=0) es- 
sentially the correlation energy calculated in the 
approximation of Brueckner and Gell- Mann.° 

We add a few remarks concerning some of the 
features and the range of application of the 
method illustrated above. 

(i) In a higher approximation Eq. (8) deter- 
mining x ,° is replaced by 


la {tH ]+31n cet =0. (10) 


as 0 


This equation represents a generalization of the 
Hartree-Fock equations to include the effects of 
correlation in the approximation of Brueckner 
and Gell-Mann.° 

(ii) It is possible to write V as a Hermitean 
form in the variables bi = 147j- Writing V in this 
form and proceeding in the same manner as 
above, one obtains results formally similar to 
those given above but representing a different 
method of approximation in the calculation of Z. 

(iii) In the approximation described in (ii) the 
equation equivalent to (8) no longer represents 


———$—_ 


the Hartree-Fock equations, but in the Fermi- 
Dirac case constitutes a set of equations which 
are generalizations of (and include) the Bardeen, 
Cooper, and Schrieffer equations of supercon- 
ductivity.* 

(iv) In the approximation described in (7) the 
term equivalent to the last term of (7) represents 
(at T=0) the correlation energy calculated in the 
approximation of the Brueckner nuclear matter 
theory.°® 

(v) In the approximation described in (ii), the 
equation equivalent to (10) determines at T=0a 
self-consistent field dependent upon the reaction 
matrix, of the type introduced in the Brueckner 
nuclear-matter theory.°® 

(vi) In some cases the solutions of (8) are de- 
generate in the sense that there is a continuum 
of solutions x,,° all yielding the same value for 
L{x,,°]. In this case special measures have to be 
taken in the application of the method of steepest 
descents. This situation corresponds physically 
to the presence of collective motions of the type 
illustrated by the rotational collective motion of 
nuclei. 

(vii) In some cases there may be two or more 
solutions of (8) representing local minima of L. 
In evaluating (7) one must of course use the so- 
lution x,,° which makes L an absolute minimum. 
However, it may be that as § is varied the local 
minima change their relative positions and the 
absolute minimum switches at some temperature 
from one to another local minimum. This pro- 
cess presumably gives a description of some 
sorts of phase transition. 

These remarks will be enlarged upon in a 
series of papers now being prepared for publica- 
tion. 








* 

This work has been supported in part by the Na- 
tional Science Foundation. 

Ton leave from the United Kingdom Atomic Energy 
Authority Research Group, Harwell. 

'R. L. Stratonovich, Doklatly Akad. Nauk S.S.S.R. 
115, 1097 (1957) (translation: Soviet Phys. Doklady 2, 
416 (1958)]. 

*R. P. Feynman, Phys. Rev. 84, 108 (1951). 

3K. A. Brueckner and M. Gell-Mann, Phys. Rev. 
106, 364 (1957). 

Bardeen, Cooper, and Schrieffer, Phys. Rev. 108, 
1175 (1957). 

°K. A. Brueckner, Phys. Rev. 96, 908 (1954); 97, 

1353 (1955). 
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SHORT -TIME MEASUREMENT OF TIME DILATION IN AN EARTH SATELLITE* 


R. S. Badessa, R. L. Kent, and J. C. Nowell 
Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received May 25, 1959) 


The purpose of this note is to describe a method 


for measuring the gravitational red shift postulated 


by Einstein’s general theory of relativity. As 
popularly envisioned, this effect causes clocks 

at different gravitational potentials to run at dif- 
ferent rates.’ The greatest difference obtainable 
in the vicinity of the earth, however, is only a 
few parts in 10*°. Nevertheless, measurements 
of gravitational red shift now appear to be feasi- 
ble. Atomic clocks with stability approaching 1 
part in 10%" have been developed, and earth 
satellites capable of orbiting at altitudes which 
provide a significant difference in gravitational 
potential with respect to the earth are available. 
The method to be described here is basically a 
short-time-measurement version of a proposed 
experiment?»® involving comparison of the time 
kept by a stable clock in a satellite and an identi- 
cal one on the ground. By “short-time measure- 
ment” we refer to the interval required to per- 
mit the small rate difference of the two clocks 
(or the frequency difference of two stable oscilla- 
tors) to accumulate a measurable difference in 
time (or phase). In previously suggested experi- 
ments the required measurement interval has 
been estimated to be several days or even weeks. 
In the system described in this note the interval 
may be as short as 30 seconds or less. 

A short-time measurement of gravitational 
frequency shift offers the following significant 
advantages: 

1, With elliptical orbits, frequency differences 
can be observed at various radii, to obtain the 
variation of frequency shift with gravitational 
potential. 

2. If the measurement is repeated when the 
satellite is at a particular radius (as at perigee), 
a slow frequency drift between satellite and earth 
clocks can be determined and removed. Since 
the quantity of interest is the variation of fre- 
quency difference as a function of altitude, a 
fixed set between clocks is inconsequential. 
Therefore, quartz crystal oscillators can prob- 
ably be employed. 

3. Continuous integration or counting over many 
orbital periods is not required. Skipped cycles 
or temporary failures would not invalidate the 
experimental results. 


The instantaneous frequency difference between 
an oscillator on the earth and an oscillator on a 
satellite as received on the earth is influenced 
by the variation in path length (first-order Dop- 
pler). Since this is typically of the order of 1 
part in 10°, and the gravitational shift is of the 
order of a few parts in 10°, a short-time meas- 
urement of gravitational shift utilizing one-way 
propagation would require extremely accurate 
knowledge of the path length. A technique for 
avoiding this difficulty is to transmit to the 
ground from the satellite a signal which has been 
“precorrected” in frequency by the increment 
required to cancel the first-order Doppler shift. 
For the following step-by-step description of 
such a system, a point of observation is selected 
at an earth pole (O in Fig. 1) in order to present 
a simple analysis that adequately illustrates the 
principles of measurement. 

1. A transmitter of frequency f is on the earth 
and transmits to the satellite. 

2. From the theory of relativity, it can be 
shown that the frequency of the signal received 
at the satellite, as observed on the satellite, is 


f'=f(1-Bcosa)/(1 ~ 37 429)2, (1) 


where £6 is the magnitude of the satellite velocity 
at the instant of reception, normalized in terms 
of the velocity of light; a is the angle (measured 
in the coordinate system of an observer at O) 





EARTH ELLIPTICAL 


ORBIT 
FIG. 1. Geometry of the measurement problem. 
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between the propagation path O-S and the velocity second-order Doppler shift, denotes a “slowing- | 




















vector 8; ¢ is the difference in gravitational down” of time in the satellite as a result of its 
potential (positive) between points S and O, orbital velocity. This term, which has been verj. 
normalized in terms of the square of the velocity fied experimentally, can be calculated easily, ong, 
of light. the orbit of the satellite is known. Subtraction gt | 
3. The satellite contains an oscillator of fre- the calculated second-order Doppler term from 
quency 2f (measured on the satellite). The re- the measured frequency difference yields the j 
ceived signal frequency f’ is subtracted from desired gravitational shift. If the transmitter 
2f by mixing within the satellite to yield a fre- operating frequency is 500 Mc/sec, the largest : 
quency difference expected frequency difference attributable to 
aff", (2) gravitational shift will be approximately 0.5 cps, | 
This frequency can be determined within + 0.01 
4. The mixture signal f’’ is transmitted back cps by measuring the accumulated phase differ- | : 
to the earth. (Reception and retransmission at ence in a 10-second interval to anaccuracyof | tr 
the satellite are assumed to be virtually instan- +36°, ob 
taneous.) When received, it has a frequency f’’’, The authors would like to acknowledge the | of 
where* : assistance given by V. J. Bates, J. W. Graham, | ot! 
a and Professor C. L. Searle. They would also 
Staf"C -F' +80¥ MA pases). (3) like to express their appreciation for the in- | . 
Substitution for f’’ yields valuable guidance of Professor J. R. Zacharias wh 
13-9 st0* and Professor G. G. Harvey. wi 
ads 1 -Scosa -B* + 
f rectal 1+8cosa | (4) a 
This work was supported in part by the U. S. Army + the 
(Signal Corps), the U. S. Air Force (Office of Scien- | 
By means of a power series expansion, retaining tific Research, Air Research and Development Com- 
second-order terms in § and first-order terms mand), and the U. S. Navy (Office of Naval Research), 5 
in ¢, we obtain ‘Lorentz, Einstein, Minkowski, and Weyl, The pre 
Principle of Relativity (Reprint, Dover Publications, obs 
fil’ =f(1+2¢ -8?). (5) Inc., New York, 1958): A. Einstein, pp. 37-59, the 
Note that the first-order Doppler term has been “2 Singer, Phys. Rev. 104, 1 (1956). a 
canceled.® The received signal frequency differs 3... H. Aller, Astron. J. 64, 45 (1959). soa 
from the ground transmitter frequency by ‘Equations (1) and (3) are alternative forms for - 
expressing relativistic Doppler shift. The c e of 
Af =f (29 -8*). (6) ms ye om so that cally one angle, a, pen: be and 
The two terms contributing to the frequency dif- considered. The CEERI One estimated to be sae 
ference are apparently inseparable experimentally. "Tene ani = Pharetin nea ¥ 
The first, 2¢, stems from a “speeding-up” of of 8 and ¢ at a single instant (viz., the instant of re- qua: 
time in the satellite as a result of the difference ception and retransmission at the satellite). There- 
in gravitational potential, i.e., the “gravitational fore, the fact that 8 and ¢ are changing cannot result - 
red shift.” The second term, -§?, known asa in an error in first-order Doppler cancellation. | a 
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EFFECTS OF 9.2-kMc/sec ULTRASONICS ON ELECTRON SPIN RESONANCES IN QUARTZ 


E. H. Jacobsen, N. S. Shiren, and E. B. Tucker 
General Electric Research Laboratory, Schenectady, New York 
(Received June 16, 1959) 


We have observed the reduction in amplitude of 
electron spin resonances at 9.2 KMc/sec caused 
by ultrasonic vibrations in the same frequency 
range. Measurements were made on resonances 
due to manganese impurity in quartz crystals’ 
and also on those due to radiation-induced centers 
inquartz. The paramagnetic resonance signal 
from the spins at one end of a quartz crystal was 
observed as a function of the power and frequency 
of the ultrasonic vibrations introduced at the 
other end, and as a function of the monitoring 
microwave power. 

The apparatus utilizes ultrasonic techniques, 
which have been reported recently,?»* together 
with a standard microwave electron paramagnetic 
resonance spectrometer using field modulation 
and synchronous detection. A quartz crystal, in 


| the form of a cylindrical rod 0.3 cm in diameter 


and $ cm long, extends between two cavity re- 
sonators in a geometry similar to that used in a 
previous experiment.* The spin resonances are 
observed in a full-wave rectangular cavity with 
the microwave magnetic field perpendicular to 
the rod axis. Sound waves are generated at the 
other end of the rod in a tunable ultrasonic drive 
cavity of the type used by Bommel and Dransfeld® 
The vibrations then travel the length of the bar 
and interact with the spins in the spectrometer 
cavity. Experiments were carried out in the tem- 
perature range of 1.5 to 2.0°K. 

Lines due to manganese in both X-cut and Y-cut 





quartz have been investigated. The paramagnetic 
centers were distributed throughout the length of 


the bar. Two gross effects can be noted in Fig. 1. 


With constant spectrometer power, the spin-re- 
sonance signal amplitude decreases with increas- 
ing ultrasonic drive power. -The second effect of 
the ultrasonic power is the shift of the peak of 

the microwave saturation curve. Increasing the 
ultrasonic power moves the peak toward higher 
spectrometer powers. Both effects are to be ex- 
pected on the basis of simple theory with the 
utrasonic vibrations considered as equivalent to 





4second microwave field (not coherent with the 
spectrometer field). 

By varying the ultrasonic frequency while keep- 
| ing the spectrometer frequency constant we can, 
| inprinciple, determine the bandwidth for the 


interaction of the vibrations with a particular 
resonance line. A typical result for one of the 
manganese lines is given in Fig. 2. The ultra- 
sonic width differs from line to line and seems 
to be smaller than the observed microwave line- 
width for the same line. 

The narrow bandwidth observed for the inter- 
action of the vibrations with the Mn spins rules 
out thermal heating of the quartz bar by the 
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FIG. 1. Derivative of absorption as a function of the 
logarithm of spectrometer power (P) for several rela- 
tive values of ultrasonic power (U) for a line due to 
manganese in quartz. Spectrometer and ultrasonic 
frequencies are identical to within a few kc/sec. 
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FIG. 2. Amplitude of derivative of absorption as a 
function of ultrasonic frequency with constant spectrom- 
eter frequency, magnet field, ultrasonic power, and 
spectrometer power. The data are for the same line 
of Mn in quartz as the data of Fig. 1. 


ultrasonic energy or by dielectric heating in the 
drive cavity as the agent producing the effect. 
Electromagnetic feed-through is also ruled out 
as a cause of the effect by observing that the 
power in the ultrasonic drive cavity has no effect 
on the resonance signal if the quartz rod is bro- 
ken or cut at a point between the two cavities. 
The dependence of the signal on the ultrasonic 
power for manganese resonances in quartz can- 
not be expressed by a simple functional relation- 
ship, and in fact, is not the same for all the 
lines observed. This is not unexpected, since 
all the observed lines saturated with spectrom- 
eter power in a way which indicates partial in- 
homogeneous broadening® and the relationship 
between signal amplitude and ultrasonic drive 
power is affected by inhomogeneous broadening 
as well as by nonuniformity of the ultrasonic 
field.” These same factors will affect the rela- 
tion between the observed vibration- spin band- 
width (Fig. 2) and the true vibration-spin band- 
width. This true bandwidth is expected to be the 
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homogeneous width of the resonance line. 

The narrow line observed for the vibration-spin 
interaction in the Mn case is evidence against the 
spreading of the phonon modes in the immediate 
vicinity of the ultrasonic frequency. 

The effect of the microwave phonons on the 
“smoky quartz” resonance centers®;® induced by 
electron, neutron, *° or x irradiation is rather 
different. The neutron-induced centers, like the 
Mn, extended through the whole length of the bar 
whereas the electron- and x-ray-induced centers 
were generated only in that part of the bar used 
in the spectrometer cavity. As before, increas- 
ing the ultrasonic power decreases the signal 
amplitude, but in contrast with the Mn, the ob- 
served bandwidth is greater than our tuning range 
of +500 Mc/sec. 

In an attempt to understand the bandwidth of 
over 1000 Mc/sec for the vibration- spin interac- 
tion in the radiation-induced centers, we have 
carried out pure piezoelectric-acoustic experi- 
ments'»*»* on the same quartz bars. The results 
showed that the X-cut bars of the irradiated 
quartz give more echoes than X-cut bars contain- 
ing manganese. This seems to eliminate the 
possibility of a broad phonon band being respon- 
sible for the large bandwidth. It has been im- 
possible to obtain piezoelectric-acoustic meas- 
urements on Y-cut crystals because of coherence 
problems due to coupling to other modes. 

In all of the vibration-spin interaction experi- 
ments it was found that a longitudinal ultrasonic 
wave (axis of the bar parallel to x axis of the 
crystal) produces an effect similar to that caused 
by predominantly transverse waves (axis of the 
bar parallel to the quartz y axis). 

It is a pleasure to acknowledge numerous dis- 
cussions with Dr. W. Kanzig and Dr. T. G. Cast- 
ner. T. G. Kazyaka helped take most of the data. 





‘The spin Hamiltonian of the manganese in quartz 
system is not yet known. : 

2H. E. Bommel and K. Dransfeld, Phys. Rev. 
Letters 1, 234 (1958). 

*E. H. Jacobsen, Phys. Rev. Letters 2, 249 (1959). 

‘N. S. Shiren and E. B. Tucker, Phys. Rev. Letters 
2, 206 (1959). 

‘H. E. Bémmel and K. Dransfeld, Phys. Rev. Let- 
ters 2, 298 (1959). 

°T. G. Castner, Jr., Bull. Am. Phys. Soc. 4, 21 
(1959), and Phys. Rev. (to be published). 

"In the case of ultrasonic vibrations interacting with 
nuclear spins, effects attributed to standing waves have 
been reported by Jennings, Tanttila, and Kraus, Phys. 
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Rev. 109, 1059 (1958). 
8Griffiths, Owen, and Ward, Nature 173, 439 (1954). 
‘Mm. C. O’Brien, Proc. Roy. Soc. (London) A231, 404 





(1955). 
10The neutron irradiation was supervised by Dr. R. A. 
Weeks, at Oak Ridge National Laboratory. 





EXCITATION OF HYPERSONIC WAVES BY FERROMAGNETIC RESONANCE 


H. Bommel and K. Dransfeld 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received June 18, 1959) 


We have observed the emission of transverse 
hypersonic waves from a Ni film, deposited on a 
quartz surface and excited to ferromagnetic re- 
sonance. The acoustic frequency was equal to 
the ferromagnetic resonance frequency. The Ni 
film, about 18.000 A thick, was evaporated in 
vacuo on to one end face of a cylindrical single- 
crystalline quartz rod of 3-mm diameter and 12 
mm long. Both end faces of the rod were polished 
optically flat and parallel. The rod was an orien- 
ted AC-cut crystal so as to allow the transmis- 
sion and detection of transverse sound waves.’ 

As shown in Fig. 1, the quartz rod was with its 
two ends part of two cavities C, andC,. One side 
carrying the Ni film was in the cavity C, exposed 
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FIG. 1. Experimental arrangement for detection of 
hypersonic waves, generated by a Ni film in ferro- 
Magnetic resonance. 


to an rf magnetic field parallel tp the Ni surface, 
while the other side of the rod was in the electric 
field of cavity C,. Both cavities were tuned to 
the same frequency of 1000 Mc/sec. C, was 
coupled to a pulsed rf generator and C, to a re- 
ceiver and scope. A variable dc magnetic field 
was applied parallel to the rod axis. 

Acoustic waves emitted from the Ni film into 
the quartz rod were received in cavity C,. The 
received power is plotted vs the magnetic field 
in Fig. 2. The maximum acoustic power of about 
1 mw was emitted at a field at which ferromag- 
netic resonance was observed simultaneously in 
cavity C,. A small subsidiary line is visible at 
fields below the saturation magnetization of 6100 
gauss. The acoustic character of the received 
signal was evident from the delay time between 
the power pulse sent into cavity C, and the re- 
ceived signal from C,, corresponding to the 
acoustic delay in the quartz rod. Since the rod 
has AC orientation, the received acoustic waves 
were of transverse polarization. 

We have also observed that, with the same 
efficiency, ferromagnetic resonance can be 
excited by incident hypersonic waves: A signal 
of the same magnitude could be seen, if the roles 
of the two cavities were exchanged, putting now 
power into cavity C, and receiving from C,. In 
this case transverse hypersonic waves were ex- 
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FIG. 2. Magnetic field dependence of hypersonic 
power emitted from Ni film. 















VoLUME 3, NUMBER 2 





PHYSICAL REVIEW LETTERS 


JULY 15, 1959 








cited in C, and the Ni film served as an acoustic 
receiver. In both experiments the orientation of 
the de magnetic field for maximum signal was 
found to be critical within one degree. 

A uniform precession of the magnetization in 
the film is related, by magnetostriction, to a 
stress ellipsoid also precessing around the dc 
magnetic field. If the direction of the dc field 
is normal to the film it can be shown that, at 
ferromagnetic resonance, the surface will under - 
go a rotational shear motion around the dc field, 
and the film will therefore emit circularly polar- 
ized transverse sound waves into the quartz, 
on which it is plated. Since the film, used here, 
was not much thicker than half an acoustic wave- 
length, the shear strain is approximately uni- 
form over the sample. Assuming further that 
the film is acoustically matched to the quartz, 
the acoustic power emitted from the film should 
be of the order 

P=pV,2'¢?, (1) 
where p=density of Ni, V;,°=transverse sound 
velocity in Ni, A =magnetostriction constant, and 
% =precession angle. According to this estimate, 
the observed maximum power output of about 1 
mw could be accounted for by a precession angle 
of 0.2°, under the assumption that the magneto- 
striction A is at these frequencies the same as 
in static experiments. It may be mentioned, that 
for other orientations of the dc magnetic field it 
should also be possible to excite sound waves of 
twice the frequency of the ferromagnetic reso- 
nance. 

It is interesting to consider also the case that 
the film is thicker than the acoustic wavelength, 
but not exceeding the penetration depth, by using 
for example a thicker ferrite disk. When such 
a sample is excited to a uniform spin precession, 
the corresponding magnetostrictive stress will 
also be uniform and a displacement can, initially, 
only take place at the surface. Similarly to the 
piezoelectric case, elastic waves are therefore 
excited at the surface. For a uniform precession 
(precession angle ¢) acoustic power is emitted 
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from the surface according to (1), if standing 
waves and reflections are neglected. 

If the spin precession is not uniform, acoustic 
waves can also be excited by spin waves as 
treated, for example, by Kaganov et al. and by 
Kittel.? Although the spin wave excitation is 
quite different from the surface excitation for 
large samples, there is little distinction between , 
both processes if the sample is thinner than an 
acoustic wavelength. Kittel? suggested that it 
may be possible to generate acoustic waves by 
using thin ferrite crystals. The magneto-acoustic 
resonance in yttrium iron garnet, observed by 
Spencer and LeCraw,° occurred at frequencies 
orders of magnitude below the ferromagnetic 
resonance frequency and is therefore a different 
process. These experiments also demonstrate 
a new method for the excitation and detection of 
hypersonic waves, and for their transmission 
into other solids without bonding them to piezo- 
electric crystals. The efficiency of this method 
is, for Ni and at room temperature, not smaller 
than in the piezoelectric excitation. There should 
be no difficulty in extending these experiments 
to higher frequencies or other ferromagnetic 
materials. For example, we have also observed 
at 1000 Mc/sec the generation of sound with a 
much thicker disk of yttrium iron garnet bonded 
to an AC-cut quartz crystal, between room ten- 
perature and 20°K. 

We wish to thank J. F. Dillon and E. G. Spencer 
for interesting discussions and for supplying us 
with the ferrite samples, A. Grieco and P. 
Andreatch for making the acoustic bonds, E. M. 
Kelley for the evaporation of the films and, 
particularly, A. DiGiovanni for his great help 
with all the experiments. 
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234 (1958); 2, 298 (1959). 
*N. I. Kaganov and Z. M. Tsukernik, Zhur. Eksptl. 
i Teoret. Fiz. 36, 224 (1959); C. Kittel, Phys. Rey. 
110, 836 (1958). 
. G. Spencer and R. C. LeCraw, Phys. Rev. 
Letters 1, 241 (1958). 



















959 


tic 


ller 
ould 





ved 


led 
2m- 


ncer 





VoLUME 3, NUMBER 2 


PHYSICAL REVIEW LETTERS 


Jury 15, 1959 








NUCLEAR SPIN RELAXATION TIME IN SUPERCONDUCTING ALUMINUM 


Alfred G. Redfield 
International Business Machines Watson Laboratory at Columbia University, New York, New York 
(Received June 5, 1959) 


The nuclear spin relaxation time 7,, in super- 
conductors was first measured by Hebel and 
Slichter’»? just below the transition temperature 
in aluminum, and by Reif’ well below the transi- 
tion temperature in mercury. These measure- 
ments indicated that T,, decreases rapidly just 
below T,, and then rapidly increases again well 
below the transition temperature, becoming 
greater than the normal state relaxation time at 
very low temperatures. Hebel and Slichter pointed 
out! that the decrease in 7, near To. could be eas- 
ily explained only in terms of an energy gap model 
of superconductivity. They supported their ex- 
planation with a detailed calculation’ based on 
the Bardeen, Cooper, and Schrieffer (BCS) theory 
of superconductivity. Comparison of the temper- 
ature dependence of 7, as contrasted with that of 
acoustic attenuation provides confirmation of the 
spin-dependent correlation of the wave functions 
used in the BCS theory.? 

Preliminary measurements of 7, in aluminum 
by Anderson and the writer*»* were in reasonable 
agreement with the Hebel-Slichter-BCS theory; 
accuracy of these measurements has now been 
improved using a He® cryostat. The experimental 
method was basically the same as that used pre- 
viously. The sample was soaked in a 3000- 
gauss field for a time long compared to T,, so 
that the spin magnetization built up to its equi- 
librium value in this field. The field was then 
turned off in about 1/10 sec and left off for a 
variable time 7; during this time the sample is 
in the superconducting phase and the spin polari- 
zation relaxes toward zero. The field is then 
turned on to 1000 gauss and within 2/10 second 
the spin polarization is measured in the normal 
state by applying suitable rf and sweep fields. 

The sample was immersed in liquid He* to 
eliminate temperature changes due to the mag- 
neto-caloric effect; temperature was inferred 
from measurement of the magnetic susceptibility 
of a solid cylinder of iron ammonium alum, with 
an estimated error of 2%. 

The recent data are plotted in Fig. 1, together 
with earlier data® taken in this laboratory, as 
a function of T_/T (assuming T,.=1.178°K). The 
signal to noise ratio was 10 or better, and the 
signal appeared to decrease exponentially with 
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FIG. 1. Experimental and theoretical nuclear spin 


relaxation time in pure aluminum plotted as a function 
of inverse reduced temperature. 


7 for all the points shown. We do not believe that 
possible impurities affected the data of Fig. 1, 
and we are currently attempting to confirm this 
by studying the effect of known impurity concen- 
trations. 

A possible source of error in this experiment 
is the presence of regions of trapped flux. Such 
regions would give a contribution to the observed 
T, characteristic of the normal state at H=H_, 
which might not be resolved with the poor signal 
to noise ratio in the present experiment. We have 
plotted in Fig. 1 the semiempirical temperature 
dependence of 7, in the normal state at H,, based 
on actual measurements at 0.5° and on measure- 
ments of 7, above T,, as a function of field.® 
Since there is no noticeable correlation between 
the variation of the normal state 7, and the super- 
conducting T,, and since they are very different 
over much of the range, we believe that the effect 
of trapped flux is small. 

The data of Fig. 1 are in reasonable agreement 
with the earliest measurements! of Hebel and 
Slichter and of Reif, and are in detailed agree- 
ment with the predictions of Hebel and Slichter 
based on the BCS theory. For comparison we 
have plotted in Fig. 1 recent calculations of 7,, 
by Hebel,® based on the expressions found in 
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reference 2 but evaluated more accurately. 
There are two unknown parameters in the theory: 
a multiplicative constant which is determined by 
the requirement that 7, be continuous at the 
critical temperature, and a parameter 7=e,(0)/A, 
where €,(0) is the half gap width of BCS at 0°K 
(3.5k7,,) and 4 can be regarded as the half-width 
of the electronic energy levels. In reference 2, 
A is introduced into the theory by smearing or 
folding the BCS density of states function with a 
square function of width 24 and height (2A)~*. 

A can be interpreted either as an uncertainty 
width due to electron scattering mechanisms not 
accounted for in the BCS theory, or as a smear- 
ing of the BCS density of states due to anisotropy 
of the energy gap.’ 

The “uncorrected” curve in Fig. 1 is that cal- 
culated by Hebel directly from the BCS theory 
for r=15; for other values of 7 the curve is 
nearly parallel to this one for T,,/T>1.5. The 
“corrected” curves are calculated using the 
energy gap found experimentally by Biondi and 
Garfunkel from microwave measurements? [i.e., 
€,(T) proportional to, but 7.5% smaller than, the 
BCS prediction]. According to Hebel,® the cal- 
culation of T,, can be corrected for an energy 
gap different from the BCS prediction if the form 
of the BCS wave functions is correct, and only 
their energies are assumed to be incorrect. In 
that case T,,7 is a function of €,(7)/kT alone. 

Agreement with theory is excellent assuming 
r=15 [A=e,(0)/15]. This agreement is in con- 
trast with the Knight-shift observations of Reif 
and of Androes and Knight. Our measurements 
are made on particles large compared to the 
coherence length (10™° cm). The parameter A 
is nearly temperature independent, consistent 


i, 


with the idea that it represents either the inverg 
electronic lifetime against boundary scattering, 

or more likely anisotropy of the energy gap.” }t 

is probable that if A represented phonon scatter. 
ing, it would be strongly temperature dependent, 
contrary to observation. 

It is a pleasure to acknowledge the contributig, 
to this research of Dr. A. G. Anderson, who 
collaborated on the early measurements. The 
writer is also indebted to many other members 
of this laboratory, especially Mr. W. V. Kiseley- 
sky and Dr. H. A. Reich. Finally he wishes to 
thank Dr. L. C. Hebel and Professor C. P. 
Slichter for several interesting discussions and 
for communicating the progress and results of 
their research prior to publication. 
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prior to publication. 
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OPTICAL MASER DESIGN 


J. H. Sanders* 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received June 2, 1959) 


Schawlow and Townes’ have discussed the pos- 
sibility of operating a maser in the infrared or 
visible wavelength region. They point out that 
the difficulty of maintaining an adequate population 
excess in the upper of the two levels between 
which maser action takes place increases rapidly 
with the operating frequency because of the large 
rate of spontaneous emission. They envisage only 
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the use of photon excitation for populating the ap- 
propriate level in the working medium. 

In observing the resonance absorption of a 
spectral line, there exists the well-known diffi- 
culty of obtaining a large amount of light froma 
source without broadening the line to the point 
where only a fraction of the total is absorbed. I 
view of this, it appears difficult to employ effi- 
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ciently optical excitation of a maser working in 
the visible region. The use of electron impact 
excitation of the working medium would reduce 
the practical difficulty of obtaining a sufficient 
density of excited atoms, and the use of a dis- 
charge in the medium itself is suggested. Ina 
medium containing a single gas or vapor, the 
possibility of obtaining the required excess pop- 
ulation in the upper of two levels exists in cases 
where there is a favorable lifetime ratio, i.e., a 
long-lived upper state and a short-lived lower 
state. As an example of this, the 6678A line in 
helium I is due to a transition between the 3'D, 
state, mean life 1.5x 10° sec, and the 2'P, 

state, mean life 4.3 x 10-” sec. Many processes 
exist in a discharge which may disturb the relative 
population of such levels from the value implied 

by these lifetime values, and it may be possible 

to choose conditions which maintain a high pop- 
ulation in the upper state. 

Detection of the onset of continuous oscillation 
in the optical maser working in the visible region 
would be particularly simple. The Fabry-Perot 
étalon, used as the “cavity” of the maser, gives 
a concentric fringe system at infinity when the 
light source is placed between the plates, similar 


to the conventional case when the étalon is used by 
transmission. Using highly reflecting plates, of 
which only one need transmit a small percentage 
of the light, and a spacing of several centimeters, 
the width of the observed rings would normally be 
larger than the instrumental width. If maser os- 
cillation were taking place, due to the presence 
of an adequate excess population in an upper level 
in the medium between the plates, the light due 

to this process would have an extremely small 
frequency width, and lines of the instrumental 
width would appear superimposed on the broader 
lines due to spontaneous emission. To confine 
oscillation to only one order, or mode, the ap- 
erture of the étalon could be restricted to reduce 
the number of effective beams in those orders in 
which the light is not normal to the plates. The 
order in which maser oscillation would be most 
likely to occur would then be the central, or 
normal, order. 





“On leave from the Clarendon Laboratory, Oxford, 
England. 

1A, L. Schawlow and C. H. Townes, Phys. Rev. 112, 
1940 (1958). 





POSSIBILITY OF PRODUCTION OF NEGATIVE TEMPERATURE IN GAS DISCHARGES 


A. Javan 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received June 3, 1959) 


In a recent paper’ Schawlow and Townes have 
discussed a possibility for obtaining maser action 
in the optical region. In their proposed scheme, 
negative temperature is obtained by optical pump- 
ing. One may expect also that under favorable 
conditions the excitation of atomic levels by elec- 
trons in a discharge can lead, in principle, toa 
state of negative temperature. However, severe 
restrictions exist if densities of the excited atoms 
as large as those needed for maser action are 
required. The present Letter considers briefly 
these limitations and certain types of systems 
which appear to be most favorable for practical 
application of this proposal. Pure gases behave 
quite differently than certain kinds of gas mix- 
tures. First let us consider the former case. 

For the purposes of rough estimates of various 
discharge conditions as described below, let us 
make a simplifying assumption that the main 


source of population of an excited state is due to 
collisions of the first kind between the electrons 
and atoms in the ground state leading directly 
to the excited level under consideration. It can 
be shown that, at least in cases discussed below, 
other details such as cascade processes and 
collision of electrons with other excited atoms 
do not appreciably effect our order -of -magnitude 
estimates. 

In the following, the indices 0, 1, and 2 refer 
to the ground, the lower, and the upper states 
of the maser levels, respectively. If 1 /90; is 
the rate of excitation of the ith excited state by 
electron impact with the atoms in the ground 
state and 1/859 is the rate for the reverse process, 
for a Maxwellian distribution of electrons at a 
temperature T,, the relationship 


904/80 = exp(E ./kT .) 
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will hold;? where E; is the energy of the ith ex- 
cited state. For @,, much shorter than the 
radiative lifetimes of the ith excited state, T;, 
it can be shown that N;/No =exp(-E; /kT,), which 
means a Boltzmann equilibrium for population 
distribution of the excited state at the electron 
temperature T e For 6;9>T;, however, one 
obtains 


N,/No = (7/6 ,,)exp(-E ,/kT ,). (1) 
The ratio T;/% 59 is a measure of departure from 
the Boltzmann equilibrium. Such a departure is 
necessary for the problem under consideration. 
The factor 1/6 0 is proportional to the electron 
density in the discharge. It is therefore clear 
that excessively large electron densities should 
be avoided in order to satisfy 639% T;, S80 that 


N2/N, =(T2/T,)( 90/629) €xP[-(E2 -E,)/kT,) 


may be larger than unity, giving a negative tem- 
perature. 

If the nature of the excited levels to be used for 
the maser action permits, it is very desirable to 
operate at a large gas pressure consistent with 
practical considerations. This is evident since 
for a given number of atoms in an excited state 
an effective temperature as determined from the 
ratio of this number to that of the atoms in the 
ground state is lower at larger pressure. 

At a large pressure, however, the trapping of 
resonant radiation becomes an important factor. 
Consider a level which is optically connected to 
the ground state. At a pressure where the reso- 
nant photons arising from this transition are 
trapped, the effective radiative lifetime of this 
level will be increased and determined by its 
decay into the other available low-lying levels. 
This effect can best be taken to the advantage if 
the upper maser level is optically allowed to the 
ground state. In such a situation, in addition to 
an increased lifetime of the upper maser level 
at a large pressure, one expects also a more 
favorable ratio of 6,,/62,. This is evident since 
a level which is optically connected to the ground 
level, as a rule, has a large excitation cross 
section by electron impact with the atoms in the 
ground state. Furthermore, due to the nature of 
the electrical dipole matrix elements, only one 
of the maser levels may be optically connected 
to the ground level. 

For a typical example of levels which allow 
operation at a large pressure, let us consider 
the case of Ne. The 2s, and 2p,, (in Paschen nota- 
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tion)® considered as the upper and the lower mage, 
levels, respectively, meet the required condi- 
tions. At a pressure where the ultraviolet photon; 
arising from the transition 2s, to the ground state 
are completely trapped, the lifetime of this state 
is estimated to be at least a factor of 60 larger 
than that of the 2p,, level. The latter level is 
known to have a mean life of 10~ sec due to its 
radiative decay into some of the lower excited 
states. Let us consider 10 mm of Hg as a typical 
gas pressure. For production of N, =10*°/cm’, 
one obtains from Eq. (1), T,/@.,=3x10™, assum- 
ing an electron temperature at kT, =2.9 ev. The 
cross section for electron excitation of the above 
Ne level is not known. Let us assume a reason- 
able value of 0=10~*” cm? for an average of this 
cross section over the electron energy distribu- 
tion. From the above value of T,/6,, and 
T,=6x10~ sec, one obtains an electron density 
of ng =5 x10" /em’. 

The numbers obtained above are intended to 
represent only rough estimates and indicate that, 
in practice, the necessary discharge conditions 
are obtainable. Considerations such as departure 
from Maxwellian energy distribution of the elec- 
trons at the high- and low-energy limits may show 
appreciable effects. However, with view to the 
possibility of operating at somewhat of a higher 
electron temperature, no serious consequences 
are expected. 

In the above scheme, the presence of a large 
density of the metastable Ne should be avoided 
to prevent the trapping of the red Ne line which 
may lead to an extended lifetime for the 2),, 
level. Furthermore, the collision of the first 
kind of the electrons with the Ne metastables 
may give rise to an unfavorable population of the 
excited levels. This can be prevented by intro- 
ducing a very small amount of a quenching gas 
such as argon to decrease the lifetimes of the 
longer lived excited atoms such as the metasta- 
bles. 

It may be of interest to compare the above sys- 
tem with a case which does not permit operation 
at large pressure. In He the 3'D and 2'P levels 
have a lifetime ratio of T,/T, =35. The lower 
level, 2'P, has a short lifetime due to its radi- 
ative decay into the ground level. An estimate 
of 6,,/62) indicates that this ratio may be as large 
as 1/15. Thus, a negative temperature may still 
be expected if appreciable excitation of these 
levels by electron impact with the metastable He 
can be prevented. On the basis of the known 
cross section for electron excitation of the 3'D 
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level and manipulations similar to those pre- 
sented above, the following rough estimate of the 
discharge condition may be obtained: At a pres- 
sure of 5x10~* mm and for kT, =5 ev, the elec- 
tron density for production of N, =10°/cm* is 
roughly n,, =10** or about 1% ionization of the 
total number of atoms. Althoygh.such a condi- 
tion is not outside of the range of practical possi- 
bility, however, it indicates the difficulties 
encountered in this case. It is likely that a more 
favorable situation of this type may be found in 

a different atom with a larger electron excitation 
cross section, thus resulting in a reduced num- 
ber for the required electron density. 

An alternative scheme is the use of transfer of 
excitation between excited states of two different 
atoms in a gas mixture. Consider a long-lived 
state of an atom (such as a metastable state). 

This state can be populated appreciably at mod- 
erate electron densities. If an excited state of 

a second atom happens to lie very close in energy 
to that of the level of the first atom, a large cross 
section is expected to exist for an inelastic col- 
lision resulting in a transfer of excitation from 
the metastable state to the excited state of the 
other atom and vice versa. Due to the nonadia- 
batic nature of the process of collision, the levels 
of the second atom which differ in energy con- 
siderably from that of the metastable level of the 
first atom do not show appreciable cross sections 
for transfer of excitation. 

If the rate at which such transfer of excitation 
takes place is larger than the rate of the radia- 
tive decay of the excited states, the metastable 
level serves to increase the contact between the 
temperature of the excited level of the other 
atoms and that of the electrons. In this case, it 
can be shown that (G*N*/GN) =(g*n*/gn) f, where 
the capital letters refer to.one atomic species 
and the small letters to the other. The letters 
with asterisks denote the excited states, and 
those without asterisks, the ground states; G and 
gare the statistical weights of the levels involved 
and N and m are the populations of the levels. The 
factor f in this equation is given by f =e4£ / kT, 
where T is the temperature of the gas as deter- 
mined by thermal kinetic energies of the atoms 
and AE is the resonance defect and is the differ - 
ence between the excitation energy of the level 
 _— to n* minus that corresponding to 


This relationship does not hold if AE is appre- 
tlably large, resulting in reduced rate of trans- 
fer of excitation. 


In order to take full advantage of the above 
effect, large differences in the partial pressures 
of the two atoms are desirable. Otherwise, due 
to reduction in electron temperature and various 
quenching mechanisms, the population of the 
metastable level will be decreased. 

As an example, a mixture of Kr and Hg with 
Kr at about 10 to 50 mm of Hg pressure and Hg 
at about 10~* mm of Hg pressure can be used for 
obtaining an abnormally large population of the 
9'P and 6'F of Hg. These two levels lie very 
close to 5s, metastable Kr. A negative temper - 
ature at the transition originating from these 
levels to 6*D may then be expected. It is expected 
that the levels 6'D and 6’P may also show a nega- 
tive temperature because of the cascade transi- 
tion from the 9'P and 6'F into 6D level. 

The °S, metastable of He also lies in energy 
fairly close to the excitation energy of the upper 
maser level of Ne discussed above. The presence 
of a partial pressure of He is expected to enhance 
considerably the negative temperature in the levels 
of Ne. 

The transfer of excitation of the type described 
above may play an appreciable role within the 
levels of the same atom. An important example 
of this is expected to occur in the levels of Ne. 
Let us consider the group of four levels 2s,, 2s,, 
2s,, and 2s,. The level 2s, is the one emphasized 
in the above for the upper of the two maser levels. 
These four levels all fall fairly close in energy. 
The level 2s, is also allowed for an optical transi- 
tion to the ground state. The transfer of excita- 
tion within these levels is expected to result in 
particular built up of large population in the 
level 2s,, this level having the lowest energy 
within this group. Thus, an even more favorable 
transition in Ne appears to be the 2s,—2),,. This 
transition lies at 10343 wave number. 

Details of the above proposals will be published 
upon their experimental verifications. 

I would like to acknowledge helpful discussians 
with Dr. W. R. Bennett, Jr., Dr. S. J. Buchsbaum, 
Professor C. H. Townes, Dr. J. P. Gordon, and 
Dr. A. L. Schawlow. 





‘A. L. Schawlow and C. H. Townes, Phys. Rev. 
112, 1940 (1958). 

*See, for instance, Rudolf Ladenberg, Revs. Modern 
Phys. 5, 243 (1933). 

‘In the LS designations the levels 2s, and 2p; cor- 
respond to (2p°4s) *P; and (2p°3p) °S,, respectively. 
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LAMB SHIFT AND ENERGY LEVEL OF A TIGHTLY BOUND ELECTRON 


D. F. Mayers 
University Computing Laboratory, Oxford, England 


G. E. Brown and E. A. Sanderson 
Department of Mathematical Physics, University of Birmingham, Birmingham, England 
(Received June 1, 1959) 


The K-absorption edge in mercury has been 
calculated to within terms of order a?mc? (where 
a’mc? =2 ry), with the object of verifying the 
presence of the Lamb shift, which is of order 
amc*. Here, we consider Za to be of the order 
of unity, since Z is large. 

Earlier rough calculations’ to within order 
amc? of all effects other than the Lamb shift gave 
a value of 6129 ry for the K-absorption edge and 
left an apparent discrepancy of only several Ryd- 
bergs with experiment to be accounted for by the 
Lamb shift. However, the experimental values 
were incorrectly quoted in this article, and the 
actual discrepancy was ~20 ry. 

To improve on the calculation of effects of or- 
der mc? and amc? other than the Lamb shift, a 
relativistic self-consistent-field calculation with 
exchange (Hartree-Fock calculation) was carried 
out by one of us (D.F.M.). Just as in the nonrela- 
tivistic case,” the energy of the electron system 
can be expressed in terms of one-electron radial 
integrals; complete expressions have been given 
by Swirles® and Grant.‘ To within order amc’, 
the K-absorption edge is given by €;, + 5€1., 
where €1, is the eigenvalue of the 1s state, and 
5€;, is the shift arising from the magnetic in- 
teraction between the K electrons. The eigen- 
value was found to be 6152 ry, and 5¢€,,, -15 ry, 
giving a net value of 6137 ry. This is in good 
agreement with similar calculations by Cohen,° 
who obtains 6138 ry. 

Inclusion of the Lamb shift completes the ef- 
fects of order amc*. This calculation was car- 
ried out® without employing an expansion in Za, 
giving a result of -38 ry when combined with that 
of the vacuum polarization.’ The resulting theo- 
retical value of 6099 ry is to be compared with 
the experimental value® of 6107.7+0.6 ry for the 
K-absorption edge, showing that inclusion of the 
Lamb shift markedly improves the agreement. 

Corrections of order a*mc? are contributed by 
the rearrangement energy and by the finite size 
of the nucleus. The former was evaluated by re- 
peating the relativistic self-consistent -field cal- 
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culations for the ionized atom, Hg*(1s)~*, giving 
a result of 4 ry. 

The correction from the finite size of the nu- 
cleus was calculated using® 


R, 2 
G,25y | ae G,2y | R, | ‘9 (G*+F*)(5E -5V) 
R, : 


+ Go?(E+m- vitor} dr, (1) 


where F, G, and F,, G, are the perturbed and 
unperturbed Dirac radial functions, multiplied 
by r, respectively, y=F/G, y,=F,/G,, and 
dy=y-¥5- Here, R, was chosen to be a radius 
just outside the nuclear radius, and R, to be an 
appreciable fraction of the radius of the K shell. 
The 6y was calculated by numerical integration 
of the differential equation from which Eq. (1) 
was obtained. The correction to the energy level 
from a nucleus with shape represented by an in- 
verted Fermi distribution’® was -3.9 ry, effect- 
ively cancelling the rearrangement energy. 

We believe that the discrepancy of ~8 ry comes 
mainly from the inadequacy of the self-consistent- 
field method in describing correlation in angle 
of the two K electrons, which tend to stay apart 
because of their mutual repulsion. 
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ATOMIC MOTIONS IN WATER BY SCATTERING OF COLD NEUTRONS 


D. J. Hughes, H. Palevsky, W. Kley,* and E. Tunkelot 
Brookhaven National Laboratory, Upton, New York 
(Received June 26, 1959) 


The scattering of very slow, or “cold, ” neu- 
trons is a useful technique for study of the mo- 
tions of atoms in solids and liquids. In this 
method, as developed’ at Brookhaven, a beam 
of approximately 0.005-ev neutrons, produced 
by beryllium filtration, is incident on a thin sam- 
ple. Measurement of the changes in energy at 
scattering gives directly the desired information 
concerning the atomic motions in the sample. 

As the energy changes are usually large com- 
pared with the incident energy, they can be 
measured accurately. In the case of an incoher- 
ent scatterer, such as hydrogen, a further ad- 
vantage of the cold neutron technique obtains; 
essentially all motions of the atoms can be meas- 
ured because all momentum gains are possible. 
The cold-neutron method was applied to water 
inthe present experiments because of the intrin- 
sic interest of the atomic motions and because 
the inelastic scattering of neutrons in water has 
a important bearing on its function as a reactor 


The water samples, which were extremely 
thin, ranging down to 0.2 mm, were formed by 
m aluminum sample holder containing milled 


portions of the correct thickness separated by 
supporting ribs. The neutrons scattered at 90° 
were passed through a slow chopper and the en- 
ergies measured by time of flight to detectors 
five meters distant, as already described.’ 
Careful background runs were made with empty 
sample holders and the spectrum shape of the 
incident neutrons was determined by scattering 
from a thin sheet of vanadium, which is essen- 
tially an elastic incoherent scatterer. The spec- 
trum of scattered neutrons is shown on a time- 
of-flight scale in Fig. 1, corrected for instru- 
mental effects. The results were somewhat un- 
expected in that they show a number of definite 
energy levels, rather than a continuum of energy 
changes that would be expected from a classical 
liquid. Because of the quantum nature of water 
thus revealed it appears that it is not a suitable 
material for simple application of the theory of 
scattering of neutrons by liquids developed by 
Vineyard,” on the basis of the more general 
theory of Van Hove.* 

Further investigations are being made on these 
various features, particularly as a function of 
sample temperature, and only the most signifi- 














._ ENERGY IN 10-Sev 


| FIG. 1. Neutron flux observed 
at 90° from incident beam direc- 
tion. Data are corrected for 
counter efficiency and chopper 
transmission. The arrow at 5.2 
4 x 107 gives the position of the 
Bragg cutoff for the Be-filtered 
incident spectrum. The other 
arrows give the measured neutron 
energies corresponding to hin- 
dered motions of the liquid. The 
quantum energy of these states is 
the difference between the indi- 
cated energy and 5.21073 ev, 

| the incident energy. 
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cant will be considered here. The prominent 
energy gain of magnitude 0.061 ev has been 
found much earlier by infrared measurements‘ 
and interpreted as a hindered rotation, that is, 
a rotatory oscillation of one water molecule in 
the field of its neighbors. Infrared measure- 
ments are not possible at energies much lower 
than that of the hindered rotations, whereas the 
present technique extends about a hundredfold 
lower, to the region of microwave absorption. 
At energies below that of the hindered rotation 
there are a number of peaks, which are pro- 
bably to be associated with various types of 
translatory motions of water molecules in the 
potential of their neighbors.* 

A most surprising effect observed is the ex- 
tremely small energy changes shown by the two 
peaks near the energy of the incident neutrons 
(Figs. 1 and 2), representing the transfer of en- 
ergy from and to the liquid of the same amount, 
0.7x10-* ev. This energy is much less than that 
of the other energy changes found in water and 
there is no known transition in liquid water that 
could account for it. The observed energy ac- 
tually corresponds closely to a line recently ob- 
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The details of the scattered spectrum near 
” The dashed curve shows the posi- 


FIG. 2. 
the “elastic peak. 


tion of the incident spectrum as would be measured with 


a thin vanadium scatterer. The water data are seen to 


coincide with the incident spectrum in the region of the 
sharp rise, indicating little diffusive scattering. The 
two small peaks at energies lower and higher than the 
sharp edge correspond to neutron energy loss and gain 
of 0.7x1078 ev. 
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served’ in the microwave spectrum of water 
vapor, although it would be very surprising for 
this energy transition to appear in the liquid. 

The low incident energy in the present experi- 
ment should be particularly suitable for obsery- 
ing the small energy changes arising from dif- 
fusive motions of the atoms. Vineyard has con- 
sidered the slight energy broadening of the inci- 
dent spectrum and in one simple model the en- 
ergy spread, Ae, is given directly by the diffu- 
sion constant 


Ae =2hDk,’, (1) 


where hk, represents the momentum transferred 
to the scattered neutron. The surprising finding 
of the present work is, however, that no spread- 
ing of the incident beam is observed (Fig. 2), 
although for the experimental conditions one 
would expect from Eq. (1) a spreading of 0.9x10° 
ev, which would be easily observed. A careful 
analysis of the present results shows that a dif- 
fusion spreading of 0.3 x10~° ev could have been 
detected, and hence the observed spread is cer- 
tainly less than one-third of the calculated value. 
In similar measurements with water at 45°C, 
the upper limit is even lower, about one-sixth. 

The present results indicate that the so-called 
“jump time” for diffusion, instead of being equal 
to the reciprocal of the Debye frequency, which 
is about 10-** second in water, is at least six 
times longer. This finding indicates that the 
water molecule remains in its initial position for 
a number of vibrations before jumping to a new 
position. The magnitude of the jump would then 
have to extend over several molecules in order 
to obtain agreement with the known diffusion 
constant. 

In order to clarify further the results just dis- 
cussed, additional measurements are being made 
at higher temperatures and for vapor and ice as 
well as for liquid D,O. As far as the application 
to moderation of neutrons in reactors is con- 
cerned, the small energy changes here presented 
are of little significance relative to the extremely 
prominent interchange of energy with the hind- 
ered rotation. The results also show that the 
motions in water are not represented well by a 
gas consisting of points of mass 18, an approxi- 
mation often used in reactor calculations. 

Brockhouse® has made similar measurements 
on water, with both low- and high-energy inci- 
dent neutrons. His results with high-energy in- 
cident neutrons, as expected, do not reveal the 
details shown by the present measurements. 
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Both his results and ours with cold neutrons re- 
yeal the hindered rotation but there are marked 
differences concerning the other features. His 
measurements do not show the sharp peaks that 
we have observed, particularly those correspond- 
ing to the small energy changes. Instead of the 
latter, he observes a general broadening of the 
elastic peak, which is somewhat smaller than 
that expected from diffusive motions. This re- 
sult is in sharp contrast to ours, which shows 
no evidence of broadening related to diffusive 
motions. 





“Guest scientist from University of Freiburg, Freiburg, 
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SCATTERING OF 50- to 140-Mev PHOTONS BY PROTONS AND DEUTERONS” 


L. G. Hyman, R. Ely,? D. H. Frisch, '! and M. A. Wahlig** 
Physics Department and Laboratory for Nuclear Science, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 
(Received June 18, 1959) 


The Compton scattering of photons by protons 
has been studied in several experiments’ * below 
s-meson photothreshold, and in two experiments’: ‘ 
above. We report here further work below meson 
threshold using both protons and deuterons as 
targets. 

The techniques of identifying and measuring 
the energy of the bremsstrahlung x-rays from the 
M.L.T. synchrotron after they are scattered by 
thin targets are essentially as reported in Pugh 
etal. The many changes which have been made 
in the synchrotron itself, giving it greater sta- 
bility and intensity, have been described in re- 
cent M.L.T. progress reports.” The many changes 
in construction and performance of the photon 
telescope, energy-measuring liquid scintillator, 
and especially electronics, are being reported 
elsewhere.® 

In Fig. 1 experimental and theoretical cross 
sections for scattering of photons by hydrogen 
are plotted as a function of lab energy at a c.m. 
angle of 90°. The cross sections are given per 
unit solid angle in the c.m. system for direct com- 
parison with the Dlinois points® and those of other 
groups.*~* Our results are given as a shaded 
area within which the cross section must lie in 
order to give, after appropriate folding, a fit to 
a smooth curve through our observed data. A 
discussion of this technique, which of necessity 
assumes a fairly smooth cross section, is given 
by Pugh et al.' The statistical standard deviation 
is represented by the height of the shaded area, 
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FIG. 1. Differential cross sections for scattering 
of 50-280 Mev photons by hydrogen at 90° c.m. 


but statistical errors are thereby exaggerated 
somewhat, since the data were taken in channels 
of approximately 10 Mev, but the cross section 
at most energies is averaged by the counter re- 
sponse over an energy band of about 20 Mev. On 
the other hand, a possible 10% error in absolute 
normalization is not displayed. We see in Fig. 1 
that these new data, and the results of the other 
experiments at this angle, are all in good agree- 
ment below meson threshold and fit in smoothly 
with the Dlinois data. 
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The dashed lines in Fig. 1 are the predictions of 
previous calculations. As has been pointed out 
before, * the Mathews dispersior. theory curve’° 
is clearly too high in the 160-200 Mev region. 

The meson-theoretic prediction, Karzas et al.,** 
falls somewhat too low in the region below meson 
threshold and rises far above the data in the 200- 
Mev region. 

We calculated the solid curve using for each 
polarization the coherent addition of four ampli - 
tudes: 

(1) The Klein-Nishina amplitude for scattering 
by a point Dirac particle. Note that this is not 
the Powell’? (Dirac particle with added Pauli 
moment) amplitude’*; that is, in this energy re- 
gion we will describe the anomalous electromag- 
netic properties of the proton in terms of its ex- 
cited states. 

(2) The induced magnetic dipole resonance am- 
plitude for scattering by the T7=3/2, J=3/2* first 
excited state of the proton, following Austern, * 
Feld,*® and Yamaguchi,’® Phase shifts taken from 
the Chew-Low-Chiu-Lomon"’ analysis of the 7- 
nucleon interactions were used in the single-level 
approximation.” 

(3) The induced electric dipole resonance am- 
plitude for scattering by the T=1/2, J=3/2" 
second excited state of the proton.'*"" This crude 
but honest single-level calculation was patterned 
after Feld’s “atomic model” calculation of (2) 
above. 

(4) The amplitude shown by Low” to result from 
the special coupling of a pair of photons to a pro- 
ton through the decay of a virtual 7° meson, which 
thus is a function of the (unknown) 7° mean life Tzo. 
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The sum of amplitudes (1), (2), and (4) gives 
approximately the same result (not shown) as the 
Mathews dispersion theory, which was evaluated 
prior to knowledge of the 1/2, 3/2” resonance, 
provided that Tyo > 2x10-** sec, but does not give 
a very good fit to the data. This discrepancy is 
removed by addition of amplitude (3), resulting 


in the excellent fit given by the solid line in Fig,}, 


The y-d cross section, completely in the lab- 
oratory system, is given along with the y-p cross 
section in Fig. 2(a). Note that our detector has 
insufficient energy resolution ot make possible 
the separation of coherent elastic scattering by 
the whole deuteron from incoherent quasi-elastic 
scattering by the constituent proton or neutron. 
At 90° lab and 100 Mev one expects” roughly half 
elastic and half inelastic. The dashed curve is 
constructed from Mathews’ nucleon amplitudes"® 
in the impulse approximation, ** using a Hulthén 
wave function matched to electron scattering by 
deuterium.** Judging by Mathews’ curve in Fig.1 
we expect the calculation to give an upper limit 
to the deuteron cross section. But the observed 
deuteron cross section seems to be about 1.60.3 
times the hydrogen cross section over the whole 
energy range. Since the deuterium run was 
bracketed by two halves of the hydrogen run, 
under identical conditions, we find no way to 
explain this high ratio. 

The cross sections for H and D at 50° are given 
in Fig. 2(b). The sharp rise at energies below 
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FIG. 2. (a) Differential cross sections for scattering of 50-140 Mev photons by hydrogen 
and deuterium at 90° lab. (b) Differential cross sections for scattering of 50-140 Mev photons 


by hydrogen and deuterium at 50° lab. 
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g0 Mev is very probably not Compton scattering; 
the most probable origin of this low-energy back- 
ground will be discussed below. The cross 
sections above 80 Mev are not very accurate, and 
we can say only that (1) the hydrogen data above 
90 Mev are fitted well by Mathews’ dispersion 
curve and (2) around 80 Mev the deuterium cross 
section is appreciably greater than the hydrogen 
cross section but seems to drop off as the energy 
rises near meson threshold. Because of the com- 
plexity of the deuteron problem, especially in this 
region of momentum transfers, and with the nec- 
essity of including all four types of elementary 
amplitudes listed above, no attempt has been 

made to interpret the 50° deuterium data. 

A brief study was made of the 50° low-energy 
background, which had been observed previously’ 
and ascribed to inner bremsstrahlung accompa - 
nying pair production (radiative pair production). 
But as estimated by Feynman and Gomez, ”* the 
intensity of inner bremsstrahlung is not sufficient 
to account for most of this large background. We 
therefore measured, using a carbon target, the 
number of electrons coming to the detector from 


duce the observed spectrum of gamma rays at 

low energies by bremsstrahlung in the rather 

thick beryllium shield’ in front of the counter 
telescope. The dependence of this background on 
target thickness and machine energy is in agree- 
ment with what would be expected from the number 
of wide-angle pair electrons*® produced in the 
target plus the Mott scattering of small-angle 
electrons produced both in the target and upstream 
from it. 

We conclude: (1) the observed scattering of 
photons by protons at 90° c.m. over the energy 
range 0-300 Mev is described well by a model of 
the proton as a Dirac particle plus strong pion 
interactions giving the well-known 3/2, 3/2* and 
1/2, 3/2" resonances, provided the half-life of 
the 7°-meson is greater than 10-!® second. 

(2) The scattering of photons by deuterons below 
meson threshold is more than expected. 

We thank D. Luckey, L. Osborne, and the M.I.T. 
§ynchrotron staff for generous help with the ex- 
riment; B. Feld, M. Friedman, F. Low, 

(. Nielsen, F. Villars, and V. Weisskopf for 
stimulating discussions; and the U.S. Atomic 

Energy Commission, the Office of Naval Research, 
he Office of Scientific Research, the National 
%ience Foundation, the General Electric Company, 
id the Sloan Foundation for their generous sup- 
port. 
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A FORMAL OPTICAL MODEL 


J.S. Bell and E. J. Squires 
Atomic Energy Research Establishment, Harwell, England 
(Received June 15, 1959) 


A construction is given of an equivalent one- 
body potential for the elastic scattering of a par- 
ticle incident on a complex target. We consider 
explicitly the case that incident and target parti- 
cles are identical nonrelativistic fermions and 
allow fully for antisymmetry. The method is 
inspired by that of Frantz and Mills’ and removes 
some defects therefrom; it also removes the 
basis for their proposed change in phenomeno- 
logical optical model analysis. Center-of-mass 
motion is still ignored. 

Denoting by |@) the scattering state and by 
\0) the target ground state, we define the model 
wave function as 


o(r, t) =(O19(F, t)1 a), (1) 


where v is the Heisenberg field operator of second 
quantization. For |a@) we take the state 


\@) -f' atte HY SNe, #10), (2) 


wich corresponds to a source of particles of 
energy E at the point r’; if y’ is sufficiently large 
only a plane wave actually reaches the target. 

In writing, with x =(r, ¢), 


(x) in “ate * Glx, x’), (3) 


=(01T(Px), d Te’ ))10), (4) 





Soe fats sat, (OIL, +H WOedw" )}10) /25, 


We make a diagrammatic analysis of this follow- 
ing Hubbard,’ using Wick’s theorem and regarding 
\0) as “vacuum.” The denominator has the effect 
simply of cancelling all diagrams in the numera- 
tor not linked to the terminal operators 7(x) and 
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we make no error by taking the time-ordered 
rather than the retarded product, for the surplus 
contribution depends on the possibility of absorb- 
ing a particle at r’ from the state |0) and so 
vanishes for large r’. 

We construct G by a perturbation theory where 
in zero order the real forces are replaced bya 
fictitious one-body potential, in general nonlocal, 


faratv'G, o)U(E, F)U(F’, 0). 
An S matrix is defined by 
< s(t, t’) = -iH"()S(t, t", 
s(t’, t’)=1, H(t) =e fof H 


where H, is the zero-order Hamiltonian and the 
total Hamiltonian is H,+H’. In terms of inter- 
action representation operators 7, the Heisen- 
berg operator can be written 

P(x) =S“*(t, 0)Y(x)S(t, 0), 


and we have the usual expression 


~iHgt 


(O1S(=, Ay(x)sit, ty" e’)S(t, -) 10) 
(01s(@,-=) 10) 


for ¢>t’, and similarly for t<t’, where |0) is the 
zero-order target ground state—assumed non- 


degenerate. Expanding above and below in powers 
of H’, we have 





G= 


dt,-+-dt 60! T{H’ (t,)-- “H(t )}10). (5) 





vl (x’); G is therefore the sum of linked diagrams 
only from the numerator. We call “improper” 4 
linked diagram which falls into two disconnected 
parts on the removal of some particle (as dis- 
tinct from interaction) line. Reasoning familiar 
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in field theory leads to the integral equation 
G(x, x’) = Gy (x, x’) ~iGy(x, x'")W(x'’, x’"")G(x'"", x"), 


(6) 
where repeated arguments are integrated over, 
G, is the zero-order value, and -iW is the sum 
of all proper linked diagrams —omitting factors 
for the terminal lines. Then from (3) with ¢ =0, 


O(F) =o 9(F) - iG, (E, F, F’)W(E,F’, FoF"), (7) 
where 


W(E, f, F’) =| d(t - try ett - ty 


(x,x’), (8) 
and likewise for G,. 
Now if u. are a complete set of wave functions 
for the potential U, with eigenvalues E,, 
unoce. %, (Fu, *(F’) occ. (Fu *(F’) 


1G,= 2 E-E +ie ~ E-E -i 
n n 





Since all the occupied states are of negative 
energy, for positive E the sign of the infinites- 
imal i¢ is unimportant in the last term. Thus 
* , 
u (eu (’) 


a =a (9) 


o(f) - o(F) ? 


This is just the integral equation for scattering 
by an added potential W, and therefore 


V(E, #, F’) =UG, ?’) + W(E, ?, F’) (10) 


is the total optical potential. The scattering 
amplitude averaged over an interval of energy 
can be obtained from V(E +ie, fF, fF’) with € finite. 
This and other details will be discussed else- 
where. 

In conclusion we compare this account with that 
of Frantz and Mills. They define an optical model 
wave function by projecting on to the zero-order 
rather than the real target ground state: 


¢’(F, t) =(01%(F, 0)1a). (11) 


In the diagrammatic analysis it is then natural 
to use time-ordered rather than Feynman dia- 
grams. Asa result, for example, the two dia- 
grams of Fig. 1 are regarded as distinct and 
the second as “proper.” The proper parts of 
diagrams are now connected by forward lines 


FIG. 1. Time-ordered diagrams. 


only (particles) and not backward lines (holes). 
In the equation replacing (9) the n summation 
then runs over unoccupied states only, and the 
optical potential is therefore of the form 


U+PW'’, 


where P is a projection operator on to unoccupied 
states. Even in the trivial case that the perturba- 
tion is only a change AU in the one-body potential, 
the “optical” potential based on U would have 

this form. In contrast, our development gives 

in this case the natural result V =U +AU. 

A more serious objection to the use of (11) as 
an optical model wave function is that it does not 
have an acceptable asymptotic form when inelas- 
tic scattering is possible. If In) denotes a state 
of excitation A, of the target, then asymptotically 


en -1 iK,’ 
(n|9(r, 0)1a@) TF @ , 
with 
K 7=K?-2MA . 
n 0 n 
If 
(Ol=DE (nl, 
n 
then ¢’ has an outgoing part 
-1 iK,’ 
Lepr es 


which is a superposition of inelastic as well as 
elastically scattered waves. 





'L. M. Frantz and R. L. Mills (to be published); 
Frantz, Mills, Newton, and Sessler, Phys. Rev. Let- 
ters 1, 340 (1958). 

J. Hubbard, Proc. Roy. Soc. (London) A240, 539 
(1957). 
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PROTON-NEUTRON INTERACTION AND THE (p,n) REACTION IN MIRROR NUCLEI* 


S. D. Bloom 
Lawrence Radiation Laboratory, University of California, Livermore, California 


N. K. Glendenning 
Lawrence Radiation Laboratory, University of California, Berkeley, California 


and 


S. A. Moszkowski 
Department of Physics, University of California, Los Angeles, California 
(Received June 10, 1959) 


Direct-interaction reactions have been treated 
by a number of authors’ from the point of view 
of obtaining information concerning spins and 
parities of excited states of nuclei. The purpose 
of this note is to point out a new use of this type 
of reaction—a study of the direct interaction 
itself. 

It should be emphasized that the direct inter- 
action referred to is the effective neutron-proton 
potential within the nucleus. Recent theoretical 
studies indicate that this effective interaction is 
essentially the long-range part (distances 
2107** cm) of the actual two-nucleon interaction.” 
The (p,) reactions on mirror nuclei are a par- 
ticularly promising source of information on this 
effective interaction, as will be developed below. 

Consider as an example, the mirror nucleus 
reaction 


C**(p, n)N** (1) 


going to the ground state.* [It can be shown that 
the (p,m) reaction is equivalent in the view adopted 
here for the two cases, doubly closed shells plus 
one neutron or minus one proton. For example, 
the ground-state reactions C'*(p, n)N** and 

N'*(p, n)O** should be twin reactions.*] In the 

j-j coupling shell model, the last neutron in C** 
lies beyond doubly closed neutron and proton sub- 
shells.® The (p,m) reaction is regarded as in- 
volving only the neutron beyond the C’* core. The 
interaction between this extra-core neutron and 
the incident proton is written 


v=V+V P’, (2) 


b 
where V,, and V, are in general spin-dependent 
central potentials and P’ is the isotopic spin ex- 
change operator. The antisymmetric wave func- 
tions for initial and final states appropriate to 
the reaction (1) are 


ld) = {I1, 2)n(@)v(1) - 12, 1)e(1)@)},—) 


If) = eI, 2)¥()(2) 12,1) 4(2)o0)}, (9 


where /|1,2) denotes a state with nucleon 1 bound 
by the C™ core, and nucleon 2 free, and 7 and » 
are isotopic spin functions corresponding to pro- 
ton and neutron states, respectively. 

The matrix element of U taken between the 
initial and final states (3) and (3’) is 
(f Ivlé) =(1,21V) 11,2)-(1,21V 12, 1). (4) 
Thus the matrix element of U separates into two 
parts, one of which is a “direct” matrix element, 
(V;), and the other an “exchange” matrix ele- 
ment, (V_). The “direct” matrix element refers 
to the charge exchange process in which the 
charge on the incident proton is transferred to 
the bound neutron. The “exchange” matrix ele- 
ment refers to the knock-out process in which 
the incident proton is captured, knocking out the 
extra-core neutron. 

We may ignore the exchange integral in (4) com- 
pared to the direct integral because of the poor 
overlap of the bound and free states wave func- 
tions on the one hand and on the other hand, the 
good overlap of the extra-core neutron wave 
function in C** with the extra-core proton wave 
function in N’*. It is emphasized that in mirror 
nuclei, this overlap argument should be particu- 
larly valid, because, aside from Coulomb dis- 
tortion, the bound-state wave functions that enter 
the matrix elements are identical. Another way 
of looking at this is to realize that the direct 
integral refers to forward scattering, while the 
exchange integral refers to back scattering. We 
know that if the scattering of two particles inter- 
acting through a potential is calculated, say in 
Born approximation, the forward scattering 
(small momentum transfer) is much larger than 
the back scattering (large momentum transfer). 

In the above approximation therefore, the (p,") 
reaction connecting ground states of mirror 
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nuclei singles out the isotopic spin exchange part, 
Vy}, of the neutron-proton interaction inside 
nuclei. Thus a comparison of the experimental 
cross section for this reaction with a detailed 
theoretical calculation of the cross section 

ought to provide the potential V, of (2). It is im- 
portant to note that no one has yet succeeded in 
calculating the correct absolute magnitude of the 
cross section for a direct interaction process. 
Levinson and Banerjee® have given the most com- 
plete treatment in their study of proton inelastic 
scattering from C™*. They found it necessary to 
use a direct interaction with a strength of more 
than twice the free nucleon-nucleon potentiai. 
These authors suggested that the increased effec- 
tive interaction may arise from a polarization 

of the nucleons in the target nucleus. This effect 
would of course be absent in the interaction (1), 

if it is correct to think of the target nucleus as 

an inert core plus one neutron in a well-defined 
state. Moreover the (p, ) reaction involves 
unambiguously a nucleon-nucleon interaction, 
whereas inelastic scattering by direct interaction 
may proceed by particle excitation in the target, 
or by excitation of a collective state, and these 


two modes are not necessarily easily distinguished.’ 


Thus the (p,n) reactions connecting the ground 
states of mirror nuclei are particularly suited to 
arather direct measurement of the effective pro- 


ton-neutron interaction in nuclei, or more spe- 
cifically the charge exchange part of the inter- 
action. 





‘This work was supported in part by the U. S. Atomic 
Energy Commission. 

‘Many references to the earlier work can be found 
in J. S. Blair and E. M. Henley, Phys. Rev. 112, 

2029 (1958). See also N. K. Glendenning, Phys. Rev. 
(to be published). 

28. A. Moszkowski, Office of Ordnance Research 
Technical Report No. 2, University of California at 
Los Angeles, April, 1959 (unpublished). 

3Recent measurements of the angular distribution of 
neutrons from this reaction are reported in S. D. 
Bloom and R. D. Albert, Bull. Am. Phys. Soc. 4, 

321 (1959); see also Albert, Bloom, and Glendenning 
(to be published). 

‘This point as well as some others briefly referred 
to throughout this Letter are developed in detail in a 
University of California Radiation Laboratory Report 
now in preparation. 

‘Evidence strongly supporting the extreme j-j coupling 
shell model for this case is given in M. K. Banerjee 
and C. A. Levinson, Ann. Phys. 2, 499 (1957). 

SC. A. Levinson and M. K. Banerjee, Ann. Phys. 2, 
471 (1957); 2, 499 (1957); 3, 67 (1958). 

"H. Ui, Progr. Theoret. Phys. (Kyoto) 18, 163 (1957); 
N. K. Glendenning, Ph.D. thesis, Indiana University, 
1958 (unpublished). 





ELECTROMAGNETIC CORRECTIONS TO THE B"*-N"? g-SPECTRUM RATIO* 


Murray Gell-Mann and S. M. Berman 
Norman Bridge Laboratory of Physics, California Institute of Technology, Pasadena, California 
(Received June 19, 1959) 


It has recently been suggested’ that a precise 
comparison of the 8 spectra of B’* and N” transi- 
tions to the ground state of C’* would provide a 
test of the nature of the vector interaction in 8 
decay. Let the spectrum of each transition, 
divided by the corresponding Fermi spectrum, 
be called S(E), where E is the total energy of the 
Bray. Then define 


R(E) = S(E, B)/S(E, N’”). (1) 


lh reference 1 it was shown that the conserved 
vector current theory of 8 decay’ predicts the 
result® 


R(E) =const(1+AE), (2) 
vhere A is determined by the width of the 15.11- 


Mev level in C™ for y transitions to the ground 
state and comes out 


A=1.33%+0.15% per Mev, (3) 


using the measurements of Hayward and Fuller,‘ 
subsequently confirmed by Garwin.® 

According to the more usual theory of 8 decay, 
in which the pion is not assigned any intrinsic 
B-decay “charge,” we may expect a formula 
similar to (2), but with a much smaller value of 
A; the reduction factor in A should be roughly 
the factor by which yp, - Ly is reduced if the pion 
current contributions to this quantity are omitted. 
(Here Up and py, are the proton and neutron mag- 
netic moments.) A reasonable guess is that 
Uy- Hy would become about one Bohr magneton 
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rather than almost five if the pion contributions 
were omitted and thus in the old theory A should 
be something like one-fifth of the value given in 
(3). The B'?-N*? experiment, if carefully per- 
formed over a range of several Mev, should thus 
be able to distinguish the two theories. 

It has been pointed out, however, by Morita® 
and Schwarzschild’ that the spectrum ratio R(E) 
is altered somewhat by electromagnetic effects 
that were ignored in reference 1. It is our pur- 
pose here to supply a new theoretical estimate of 
the correction factor that should be applied to 
Eq. (2) to take account of these effects. With 
respect to the processes considered by Morita, 
our numerical results are in agreement with his, 
although our method is somewhat different. 

Let us state the answer first. Instead of Eq. (2), 
we should use 


R(E)/f (E) = const[1 + (A+ 5A)E], (4) 
where 


J) =1+ tet) hy 89/042) 


4 
-3) 





(ky-Ro’) k,? -k,’*l 
with a signifying the fine structure constant, k, 
the end-point total energy for B’*, k,’ the end- 
point total energy in N’*, and with 

1 &£ E aeatel 


-— (E*-m*)” i E - (E?-m*)” 


+z (5) 


= In(2E/m) for E>>m. (6) 


The calculated value of 5A is -0.25% per Mev, 
to which we would like, for the sake of safety, to 
attach a theoretical “error” of 0.15% per Mev. 
This error, and, indeed, the electromagnetic 
correction itself, is considerably smaller than 
the effect being sought, and therefore the exist- 
ence of the correction does not appreciably 
weaken the B’*-N"” experiment as a test of the 
conserved current theory. 

In deriving this correction, we have considered 
the following effects, many of which are negli- 
gible: 

(a) The effect of the B’*-N mass difference 
on the spectra. Insofar as the spectra are al- 
lowed, this is taken care of by factoring out the 
Fermi spectra as in Eq. (1). However, there are 
small forbidden corrections to the spectra caused 
by retardation effects in the axial vector coup- 
ling.. These cancel between B™ and N"™ except 
for those terms that involve the end-point energy, 
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which is different for the two nuclei. The result. 
ing correction is a contribution to 5A of the form 


6A = (Ry - ho’ JOYE SOr? - SOFT. (1) 


(b) The effect of the B’*-N** mass difference 
on the inner bremsstrahlung. The probability of 
inner bremsstrahlung for a given electron energy 
depends on the end-point energy. The ratio of 
the correction factors in B™” and N’ is finite 
without an infrared cutoff and is given® by f(E) in 
Eqs. (4) and (5). 

(c) The difference between B™ and N™ wave 
functions caused by electromagnetic violations of 
charge independence. This alters ft values but 
changes the spectra only through very small 
changes in the already small retardation terms. 

(d) The effect of virtual transverse photons 
exchanged between § ray and nucleus. This is of 
order nucleon v/c compared to the Coulomb in- 
teraction between § ray and nucleus, which we go 
on to consider. 

(e) The effect of the Coulomb interaction be- 
tween § ray and nucleus on the retardation cor- 
rections to the axial vector matrix element [<. 
This is often considered as the sum of three 
terms—the “finite de Broglie wavelength effect,” 


the effect of finite nuclear size, and the Coulomb 


corrections to forbidden matrix elements. We 
treat them all together. 

Since Z is small, we can treat the Coulomb 
field in first order perturbation theory. Between 
the 8 decay and the Coulomb interaction the nu- 
cleus is in a virtual state which may be excited. 
It is difficult to evaluate the sum over all the ex- 
cited states and we therefore look at two opposite 
approximations: 

I. The excitation energies of the important 
nuclear levels are neglected by comparison with 
fic/R, where R is the nuclear radius. In other 
words, the nuclear motion is assumed slow com- 
pared to the velocity of light with which the § ray 
is traveling. Under these conditions the sum 
over nuclear levels can be done by closure. In 
physical terms, the electron sees the instanta- 
neous positions of the various protons that are 
sources of the Coulomb field. This approxima- 
tion should be practically valid. As a check, 
however, we also consider the opposite case. 

Il. The energies of all excited states are taken 
to be so large (in energy denominators) that only 
the ground states survive in the sum. Physically 
this corresponds to the electron’s seeing the nu- 
cleus as a smeared-out charge distribution. 
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What is important for the B’*-N” ratio is the 
charge distribution in C* and the average of the 
Bp” and N’? charge distributions. The latter is 
just the charge distribution of the excited state 
of C* at 15.11 Mev. For simplicity we assume 
that to be the same as the charge distribution of 
the ground state. ; 

In each case we calculate the terms of order 
ZaRE in the spectrum ratio, dropping higher 
orders in RE. The result in case I is a contribu- 
tion to 5A equal to 


6A,” = -4a{ fo} 2D A$ [olt-F, I+ $[OF-(F-F,)1F-F 1 


- §/ORE-F FF}, (8) 


where the sum is extended over all protons ex- 
cept the decaying nucleon. In case II the result 
is 


0A,” = -4a{fo} [oF ab r'{2 folr-F' | 


+ $Jor-(F-F’) |F-F" |"? 
- $/o-r(F-F’) |F-F"|~*}, (9) 


where p is the charge density in C™ in units of e. 

(f) The Coulomb effect on weak magnetism 
itself. Since the weak magnetic effect on the 
spectra comes from V-A interference, it is of 
opposite sign in B** and N’*. The first order 
Coulomb correction to it thus has the same sign 
in B’* and N*” and practically cancels. 

(g) The Coulomb effect on the relativistic axial 
vector matrix element fy,. Like the uncorrected 
Jy,, this has a negligible effect on the spectrum. 

We have estimated 5A™, 6A,, and 5A,,™ 
using the shell model with harmonic oscillator 


wave functions adjusted in radius to give the 
charge distribution of Hofstadter.* The results 
are that 6A = -0.036% per Mev, 5A;™=-0.21% 
per Mev, and 5A;;=-0.23% per Mev. Since 
case II is so close to case I and since case I is 
much closer to the truth, we take 6A =6A%+6A;™ 
= -0.25% per Mev. 
The largest error in our result probably comes 
from the use of the shell model. We can esti- 
mate how bad the shell model is by using it to 
calculate the allowed matrix element {6 for the 
B decay of B’*. The experimental matrix element 
is smaller by a factor of 0.4 (in absolute value— 
presumably the sign is right). The ratios in 
Eqs. (7)-(9) are probably calculated slightly better 
than the absolute value and we have therefore 
assigned an error of about 0.15% per Mev to 6A. 





* 

This research has been supported by the U. S. 
Atomic Energy Commission and by a grant from the 
Alfred P. Sloan Foundation. 
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2S. S. Gershtein and J. B. Zeldovich, Zhur. Eksptl. 
i Teoret. Fiz. 29, 698 (1955) (translation: Soviet 
Phys. JETP 2, 576 (1957)]; R. P. Feynman and M. 
Gell-Mann, Phys. Rev. 109, 193 (1958). 

3 Everywhere except in Eq. (6) we neglect the mass 
of the electron. 

‘E Hayward and E. G. Fuller, Phys. Rev. 106, 
991 (1957). 

5T. T. Lauritsen (private communication). 
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8S. M. Berman, thesis, California Institute of 
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PREDICTION OF DEUTERON POLARIZATION FROM d-a SCATTERING 


R. J. N. Phillips 
Atomic Energy Research Establishment, Harwell, England 
(Received June 12, 1959) 


There is current interest in possible means of 
producing and analyzing deuteron polarization. 
d-a scattering on the 1.07-Mev resonance’~® and 
the T(d, n)He* reaction* have been studied. The 
present note shows that d-a scattering can be a 
useful producer (and analyzer) of deuteron polari- 
zation at energies between resonances, where 





the weak energy dependence may be an advantage. 


We have used the phase shifts of Galonsky and 
McEllistrem,® deduced from differential cross 
sections with the help of dispersion theory, to 
calculate polarization parameters in the interval 
between the resonances at 1.07 and 4.6 Mev 
(deuteron lab energy). 

In the region of 2 Mev there is considerable 
vector polarization, increasing smoothly with 
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FIG. 1. d-a vector polarization near 2 Mev. 


energy (Fig. 1); the tensor polarization param- 
eters (see below) are everywhere less than 3% 
and are not shown. These effects come mainly 
from the “tail” of the *D, resonance (1.07 Mev) 


interfering with potential and Coulomb scattering. 


At higher energies, however, °D, and °D, scat- 
tering (resonances at 4.6 and ~5.9 Mev) become 
more important. At 3.5 Mev the vector polari- 
zation has approximately doubled, while tensor 
polarization is no longer negligible (Fig. 2). 
Above this point the scattering becomes more 
strongly energy-dependent and is dominated by 
the 4.6-Mev resonance. 

Our notation follows that of Stapp.® If k and k’ 
are initial and final relative deuteron momenta, 
we define unit vectors N, K, and Q in the direc- 
tions k xk’, k’ -k, and k’ +k, respectively. Then 
if S; (i =1, 2,3) are the components of the spin-1 
operator, the vector polarization P;=(S;) and 
tensor polarization T; ;= (3(S;S; +S; jS;) - 26, ) set 
up in the scattering of initially Pal 
deuterons have the forms 


P,=BN > (1) 


T, 47 (NN - 


Ps * 
35; ,) + K(@;@, KK) 


+1(QK+K,@,). (2) 








-1.o- 


FIG. 2. d-a polarization at 3.5 Mev. 


The analyzing power of the scattering is described 
by vector and tensor “analyzabilities” P; and Ti, 


P; —, ° (3) 


Ty*Ty ;- 210K, +K@)), (4) 
such that the cross section o for scattering polar- 
ized deuterons is 
= 3 

o=0,(1+ ip P. +37, T,). (5) 
Here P; and T;; describe the incident polariza- 
tion and d is the unpolarized cross section. For 
pure vector polarization, this is like the familiar 
scattering of spin-4 particles, except for the 
factor 3 in Eq. (5). 





'L. J. B. Goldfarb, Nuclear Phys. 7, 622 (1958). 

*L. J. B. Goldfarb and J. R. Rook, Nuclear Phys. 
(to be published). 

3L. G. Pondrom, Phys. Rev. Letters 2, 346 (1959). 

‘Galonsky, Willard, and Welton, Phys. Rev. Letters 
2, 349 (1959). 

5A, Galonsky and M. T. McEllistrem, Phys. Rev. 
98, 590 (1955). 
a P. Stapp, thesis, University of California Radia- 
tion Laboratory Report UCRL-3098, 1955 (unpublished). 
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AH HYPERFRAGMENT EMITTED FROM AN ANTIPROTON CAPTURE STAR* 


A. G. Ekspong, A. Frisk, and B. E. Ronne 
Institute of Physics, Uppsala, Sweden 
(Received June 18, 1959) 


The observation of a definite case of A° hyperon 
emission from an antiproton annihilation in a 
medium heavy nucleus is reported. The lambda 
hyperon is emitted bound in a , H* hyperfragment. 

It is known’ that antiproton annihilations, in 
general, produce a number of pions and that 
strange particles occur only rarely. From 
emulsion experiments, Chamberlain et al.° esti- 
mate that about 3.5% of all antiproton stars con- 
tain KK meson pairs. A similar result has been 
obtained in our study. When the annihilation 
takes place in a nucleus, one may expect some 
cases where hyperons are produced, either by a 
secondary process involving the reabsorption of 
aK meson or by a primary process involving 
two nucleons and the antiproton. Observations of 
two possible = hyperons emitted from antiproton 
capture stars have been published.'’* In view of 
the uncertainties attached to these earlier ob- 
servations, we regard our observation to be of 
some general interest. 

The event occurred in an emulsion stack ex- 
posed to a beam of 740-Mev/c antiprotons at the 
Bevatron in Berkeley. The number of antipro- 
tons observed in this experiment is 374. Details 
on the scanning technique and the identification 
of the antiprotons have been published elsewhere? 

A microphotograph of the antiproton star with 

















FIG. 1. Microphotograph 
of the antiproton annihila- 
tion at A and the mesonic 
decay of the , H‘ fragment 
atB, The kinetic energy of 
the hyperfragment is 54 
Mev. Path lengths are in 
millimeters. Track 10 is 
due to an 18-Mev negative 
™meson, track 11 to an 89- 
Mev He nucleus. All other 
prongs in star A represent 
protons. 


the hyperfragment is shown in Fig. 1. The par- 
ticle which is assigned to be the antiproton was 
found to be moving towards the star and to have 
protonic mass. As the kinetic energy of this 
particle at the star is only 96+6 Mev and the en- 
ergy release is at least 979 Mev, it is clear that 
the star represents the annihilation of an anti- 
proton. All the tracks connected with this event 
are completely contained in the stack. Each of 
the eleven outgoing tracks are numbered as seen 
in the figure. Number 7, which according to our 
measurements most probably is a proton of 150- 
Mev energy, interacts in flight. All other par- 
ticles come to rest. 

Number 2 is the hyperfragment. This track is 
very suitably located for accurate measurements. 
It is 3630 microns long and passes through 3 
plates before decaying at rest (at point B in the 
figure). The charge of the fragment was meas- 
ured by counting the number of 5 rays along the 
entire track length. The result Z =1 is certain. 
The mass was determined by the constant sag- 
itta scattering method to be 4.0+1.1 proton 
masses. Calibration results from measurements 
on protons were utilized for both the charge and 
the mass determination. These results indicate 
that the particle is probably a H® or H* hyper- 
fragment. 
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The two secondary tracks, one of which is a 
negative pion, are collinear within experimental 
errors. They are emitted at an angle of 22.5° 
with the emulsion plane (in unprocessed emul- 
sion). The range of the secondary 7” meson is 
38.7 mm and it gives an absorption star at rest 
(at point C in the figure). The kinetic energy of 
the pion is then 52.3+1.0 Mev, where the stand- 
ard error is computed from straggling error, 
observational error, uncertainty of emulsion 
density, and uncertainty of shrinkage factor. The 
other secondary has a range of (8.0+0.4) u. The 
approximate collinearity of the two secondary 
prongs suggests that the event represents a two- 
body decay. We assume that this is the correct 
interpretation and will later return to a discus- 
sion of this assumption. The nature of the re- 
coil particle is deduced from its observed range 
and its momentum, which is equal to the ob- 
served 7-meson momentum on the above-men- 
tioned assumption of a two-body decay. On the 
basis of momentum, we expect the range to be 
8.0+0.4 p if the recoil particle is He*, and 
11.5+0.5 yp if He*®. The measured range of 
8.0+0.4 yu is consistent with a He* recoil. This 
leads to the following interpretation of the event: 


where Q =54.621.0 Mev. 

Other interpretations of this event were con- 
sidered. In addition to the above elimination of 
a He® recoil because of range, it is further ruled 
out in a two-body decay, , H’~-He*+7°+Q, be- 
cause it leads to a negative binding energy of 
-12.5 Mev for the A hyperon. Although the ob- 
served collinearity indicates a two-body decay, 
we have also considered the only possible three- 
body decay, namely , H*~He*+n+2°+Q. This, 
too, is ruled out because it leads to a negative 
binding energy of the A hyperon (in this case 
By, <-15.7 Mev). 

The result is that the event is consistent only 
with a , H* fragment undergoing mesonic two- 
body decay. The binding energy of the A in , H* 





is then By =2.6+1.0 Mev, which is consistent 
also with other measurements of this quantity.’ 
In our calculation of the A-binding energy, we 
used the following input data: Q, =37.36 Mev, 
and proton separation energy in He‘ =19.80 Mey, 

A limited search for a neutral K meson in as- 
sociation with the event was undertaken. None 
was found in a volume such that the probability 
to find the decay of 6, meson was about 5%. 

We conclude that the observation establishes 
the emission of a A hyperon bound in a , H* frag- 
ment from an antiproton annihilation in a medium 
heavy nucleus. The binding energy of the hyperon 
in , H* is 2.6+1.0 Mev. 

We are greatly indebted to Dr. E. J. Lofgren, 
Professor G. Goldhaber, Professor E. Segré, 
Professor O. Chamberlain, and others at the 
Bevatron for arranging the exposure. We would 
like to thank Professor K. Siegbahn for his sup- 
port and encouragement of this work. The sup- 
port given by the Swedish Atomkommitté is 
gratefully acknowledged. 
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Interest has recently developed in the subject 
of the electric polarizability of the neutron.’»? 
This is a quantity which can presumably be cal- 
culated when one has available sufficient knowl- 
edge of the structure of the neutron, in particu- 
lar, its mesonic structure. Estimates of the 
polarizability based on current meson theory or 
on photomeson production data suggest a value 
lying between 4x10-* cm® and 2x10-® cm’, 
while currently available experiments only place 
an upper limit of about 10-“ cm® on its value. 

It is of some interest to remark that the mere 
existence of an anomalous (Pauli) moment for 
the neutron already implies a negative electric 
polarizability for the neutron, in a manner some- 
what analogous to the situation for the neutron- 
electron interaction.* This may be seen by sim- 
ply writing the Dirac Hamiltonian for a neutron 
with a Pauli moment Ln in a uniform static ex- 
ternal electric field E: 


H =6Mc*+ca-p +i, parE. 
Since p commutes with the Hamiltonian we can 


find the eigenvalues of the energy W for states 
in which p =0: 


Hy =(6Mc* +i parE)y= Wy. 
Then, since 
Hty = (M*c*+ u 3E*)p= Wy, 


one has immediately that the energy eigenvalues 
are 


Was [M*c*+ u 2E*} = [Mc*+ up 3E*/2Mc*+--], 





ELECTRIC POLARIZABILITY OF THE NEUTRON* 


L. L. Foldy 
Case Institute of Technology, Cleveland, Ohio 
(Received June 15, 1959) 


demonstrating an electric polarizability of 
Q=- H /Mc* = -6.2x10-“ cm’. 


While this is somewhat smaller than anticipated 
contributions from the polarizability of the meson 
cloud, it is questionable whether it should there- 
fore be dismissed lightly. It need only be re- 
membered that the Pauli moment contribution to 
the neutron-electron interaction was also ex- 
pected® to be small compared to the intrinsic in- 
teraction arising from the meson cloud, yet ex- 
periments have shown it to dominate the ob- 
served interaction for reasons still unfathomed. 

It is interesting that should the intrinsic polar- 
izability of the meson cloud indeed be anomalously 
small, the resultant polarizability of the neutron 
may be negative, a result which we believe not to 
have been anticipated. While it is easy to trace 
the negative sign of the magnetic contribution to 
the contribution of negative energy states in the 
one-particle theory, we have not been able to 
construct a rationalization in terms of classical 
or semiclassical concepts. 





*Research supported by the U. S. Atomic Energy 
Commission. 
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PROTON POLARIZATION IN MULTIPLE PION PRODUCTION AS A CONSEQUENCE 
OF THE ISOBAR MODEL™ 


L. F. Landovitz and L. Marshallt 
Brookhaven National Laboratory, Upton, New York 
(Received June 24, 1959) 


The isobar model'»? for the system p+z* is in 
reasonable accord with the gross features of the 
momentum distribution for pions emitted in pion 
production in proton-proton collisions, and also 
in photoproduction of pions from nucleons. In the 
present note we shall show some consequences 
which follow near threshold for production of the 
isobar by pions if the model is taken seriously, 
i.e., the isobar production is assumed to be the 
dominant inelastic process. 

We consider the process 1*+p—7 +N", followed 
by N*~1+N. The isobar N* is in a pure state 
T=%, J =%, positive parity. Because the initial 
system 7*+p) is in a pure T =3 state, we need not 
consider any T =} states. In Table I are shown 
the relevant angular momenta for initial and 
final configurations. From these it follows that 


Table I. Possible angular momentum and parity states. 





Orbital angular Total angular Orbital angular 





momentum of momentum momentum of 

incident state / J final state L Parity 
0 1/2 2 (-) 
1 3/2,1/2 1 (+) 
2 5/2, 3/2 2,0 (-) 





if d and higher angular momenta are not excited 
in the incident configuration of 7++p, and if the 
final configuration has low energy, then only the 
Py, and P,, angular momentum states can be ex- 
cited in the configuration N* +7. With these 
assumptions, the production process 1+p—N* +7 
may be described by only two amplitudes a, and 
a, corresponding to the states P,, and Py. 

We assume that the mass of the isobar cor- 
responds to the peak of the p resonance. Of 
course, because of the width of the resonance, 
it is possible to produce the isobar at somewhat 
lower energies. The laboratory kinetic energy 
of the incident pion able to produce an isobar of 
this mass is 380 Mev; however, in practice, in 
order to obtain significant production one would 
use pions of somewhat higher energy. The cross 
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i, 


section*® for inelastic processes in 7p scattering 
remains relatively small, ~3.6 mb, as compared 
with the total cross section, ~16 mb, at 400 Mey, 
from which one concludes that s and p waves 
alone are adequate to explain the magnitude of 
inelastic processes. Consequently it is reason- 
able to neglect d waves‘ in the initial system 
p+n*. In the final system the particles are mov- 









ing so slowly that d wave is ruled out here also. 

In the case in which the isobar decay products 

are emitted in the production plane, the spin of 
the decay proton can be polarized? only perpen- 
dicular to this plane. We calculate its magnitude 
here in terms of the center-of-mass angles 6, 
and 6,, where 6, is the angle between the momen- — 
tum vectors of the pions involved in the produc- 
tion of the isobar and where 6, is the angle be- 
tween the momentum of the pion produced by 
decay of the isobar and the incident direction. 
In what follows, we have disregarded the small 
change in 6, caused by the transformation from 
the rest system of the isobar to the barycentric 
system for production. 

The amplitude for the production process 
1*++p—N +n was written in terms of a, and a,, 
the partial amplitudes for J=3 and J=3, respec- 
tively, assuming the Z axis to be along the 
momentum vector of the incident pion. 


1 
F (8, g) =6° ye, 17a, - V38Uye, 17a. 


This expression was re-analyzed into states of 
the isobar of spin 3, and of the associated pion. 
The isobar then was assumed to decay into a 
m+N; correspondingly the isobar wave functions 
were rewritten in terms of wave functions for 
the final 1, N system. Finally the spin functions 
for the final proton were projected onto an axis 
perpendicular to the scattering plane, which is 
here taken to be the y axis. For events in which 
all particles lie in the scattering plane, azimuthal 
angle ~, =~, =0 where subscripts 1 and 2 are as 
before, we find the polarization of the decay pro- 
ton to be given by the expression 


. * 
_,2 (@,a," - a,a,") 
P(6,, 0) =iz (6, 6.) [6cosé,sin@,cos*6, 


+9sin?6,cos6,sin6, -3cosé,sin@, sin”6,] , 
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where 0(6,, 6.) is given by 


A significant consequence of this model is that 
the protons are not polarized for certain values 
of the angles 6, and 6,. Figure 1 shows the cor- 
relation of 6, and 6, for those values which give 
zero polarization. Of course the experimentally 
observed polarization will not be exactly zero at 
these angles because of the background direct 
or nonisobaric production of 7+7+N. But if the 
model has validity, then the polarization should 
be low or minimal in the regions indicated in 
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FIG. 1. Contour diagram of polarization. 










o(6,, 6) = |a, |?[3(cos*6,cos*é, + sin*6,sin?6,) + 1 + 6cosé@,cosé,siné@, sin8, | 


+ 21a, |?[1 +3cos?6,cos*0, + 15sin?6,cos*0, - 12cos0,cosé@,sind, sind, | 


+(2)3(a,a," +a3@, )[4cos*6,cos”6, - 3cos6,sin@,cos6,siné@, - 2sin?6,sin?6, 


+cos?6,sin?6, +4sin?6,cos*6,] . 





Fig. 1. Furthermore, it may be seen by inspec- 
tion of the expression for the polarization that 
its relative sign is different on opposite sides of 
the lines of zero polarization. It may be mer- 
tioned that there exist configurations of isobar 
production by 7*+p such that the directly pro- 
duced pion can be readily distinguished from the 
decay pion by its angle of emission provided that 
the angle and energy of the decay proton is meas- 
ured. 

We wish to thank R. Serber, -G. C. Wick, C. N. 
Yang, T. D. Lee, and M. Goldhaber for helpful 
conversations. 
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VARIATIONAL PRINCIPLE FOR GEOMETRODYNAMICS 
David Sharp 


i 


Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 
(Received June 23, 1959) 


The purpose of this note is to propose a varia- 
tional principle for geometrodynamics, to derive 
from it certain of the field equations, and to re- 
late the problem of obtaining the remaining field 
equations to the problem of finding a superpoten- 
tial for the metric tensor. 

Misner and Wheeler,’ following Rainich,” have 
shown that the content of general relativity and 
classical source-free electromagnetism can be 
fully expressed in terms of quantities depending 
only on the metric tensor g,,,, of a Riemannian 
space. In their formulation the field equations, 
fully equivalent to the usual ones, are 


a ap 
RAR pF yy Ry af )/4, (1) 
R “.R=0, (2) 
iv 
Ro 2 9, (3) 


ty Mp7 (4) 


where the vector a B is defined as 
AyluL V OT 
a =(-g)[Bryuvj]R R’/R_R. 5) 
P (-g)*[2a uv] fig = ( 


For some purposes it is useful to replace the 
10+1=11 equations (1) and (2) by the equivalent 
10 equations: 

RR” 
R--- > +4 
uy RR ‘Buy 
OT 


«_ ie 
2 


OT 
(RF )¥? = 0. (6) 


The identity of (6) to (1) and (2) may be seen by 
setting =v and summing. The factor -} could 
just as well have any other constant value —hence 
the quotes. 

By way of explanation it may be recalled that 
the curl condition (4) guarantees the existence of 
a scalar, 

x 


a(x) -/ a ax" Qo, (7) 
0 
the “complexion” of the electromagnetic field. If 
a=nn (n=0,1,2,...), we have a pure electric 
field; if a=inn (n=1,3,...), the field is pure 
magnetic. The complexion at a point, together 
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with the value of the Ricci curvature tensor at 
the point, completely determine the electromag- 
netic field f yy and hence also the stress- energy 


tensor, 
_ @ 4 OT 
a ee ( 


at that point. 

The algebraic equations (1)-(3) satisfied by the 
Ricci curvature express the two-way connection 
between the electromagnetic field tensor f uy and 
the Ricci tensor R y’ Thus, any electromagnetic 
field tensor f,,,, produces a Ricci curvature R 9 
that satisfies (1)-(3). Conversely, given anR 4 
that satisfies (1)-(3), we can solve for fu py» the 
solution being uniquely determined except for a. 
The curl condition (4) subsequently determines a 
up to an additive constant ap. 

It is often very helpful to have a variational 
formulation of a set of field equations. 

Such a formulation is well known’>* in the more 


familiar version of general relativity plus electn- 


magnetism, where the f uy are treated as foreign 
objects immersed in the metric space rather than 
expressed—as above in Eqs. (1)-(5)—directly in 
purely geometric terms. 

In the traditional version of general relativity 
plus electromagnetism, the coupled field equa- 
tions are 


ait a - oT 
‘., 8 FS oS . 8 Sot , (9) 


bv 
gat (10) 


* Vv 
fee a0. (11) 


The variational formulation of these equations 
has the form 


6/=0, 1= var, dt =dx°dx' dx*dx*, (12) 


with 


£ =(Lagrangian density for the gravitational field 
+ (Lagrangian density for electromag- 
netic field) (13) 


=R(-g)*+4(-g) a’ (14) 


The usual procedure to obtain the field equa- 
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tions from the above is to require (for detail, see 
references 3 and 4): 
(a) The existence of a vector 4-potential for 





fur’ 


f 


uv uly vip (15) 


— ,HV 
= lv " 
(0) fy (or %,) and g"” are treated as the quan- 
tities to be varied; 
(c) 59, and 5g” are both to vanish on the sur- 
face o bounding rT. 
We propose for geometrodynamics the Lag- 
rangian density 


¢ =R(-g)* - 3(-g)"(R oy wa 


sons | f "* | 








(16) 


The first term is, as before, the Lagrangian 


) density for the gravitational field. Likewise, the 


second term is an expression for the Lagrangian 
density E* - H* for the electromagnetic field, ex- 
pressed in purely geometrical form. This ex- 
pression can be derived from results in Misner 
and Wheeler.°® 

We see that our Lagrangian density is a func- 
tional of the metric tensor g#” and its derivatives 





up to the second order, and of the path of integra- 
tion s occurring in the definition of a(x), Eq. (7). 
We require that J = [¢dz be stationary with res- 





pect to arbitrary variations of both the path 7 and 


gh 
We first vary with respect tom. 6,J=0 results 
in the condition that the integral 


f ° 
a dx’, 
- 


(17) 


be true if the integrand a B has a vanishing curl: 


a, -@ _=0. (4) 
Bly ylB 
This is one of the equations we seek. We may 
how write £ in a slightly simpler form: 


£=R(-g)" - a ll cos2a. (18) 

With no loss of generality we can set a(x) =0 at 
the point of variation, since this amounts to 
nothing more than choosing a convenient origin of 
coordinates. Note, however, that this does not 
imply that 5a =0, although 5(cos2q) =0. 

At this point it would seem natural (1) to vary 
Eq. (18) with respect to g#”, (2) to integrate by 








parts terms containing derivatives of 5g", and 
(3) to set equal to zero the resulting coefficient 
of 5g4¥. The field equations obtained in this way 
are not satisfactory; they are not equivalent to 
the desired and expected Eq. (6). 

It seems reasonable to attribute this difficulty 
to a wrong choice of the quantities to hold fixed 
on the surface, not to a wrong choice of quanti- 
ties to be varied. 

This is not the first time that such a problem 
has arisen, as witness the variational principle: 


5 f (E* - H*)d(vol) = 6f (-g)*f , a iP git=0, 
(19) 


for electromagnetism alone in a space with pre- 
scribed metric g,,,. 

The most primitive analysis of this variational 
principle, based on the assumption that the f y 
are the quantities to be varied, and that the varia- 
tions of f,,,, vanish at the boundaries, leads to 
the incorrect result that f uv= 0 throughout the 
interior. 

The customary procedure is to add a supple- 
mentary condition, Eq. (15), and treat the y,, as 
the quantities to be varied, and hold them fixed 
on the boundary. 

By setting the coefficient of 5g,, equal to zero, 
we obtain four of Maxwell’s equations. The other 
four were assumed from the outside, Eq. (15). 

However, there is still a third approach, °~® 
due to Finzi, whereby one can obtain all eight of 
Maxwell’s equations without the addition of such 
a subsidiary condition from outside. 

Finzi begins by observing that any antisym- 
metric tensor f py can be expressed in the form 

(-g)*[uvap]e"%2"?B 


[als] 20) 





Fuv*Atuty]* 


where A,, and By are 4-vectors. We now vary 
the action integral with respect to both A,, and 
B,, Subject to the condition that they both be 
fixed on the boundary. The demand that the ac- 
tion be stationary immediately produces all eight 
of Maxwell’s equations. 

Judging by the analogy between the f pv 28 field 
variables in electromagnetism and the g4” as 
field variables in gravitation, we expect that 
there are alternative ways to vary arbitrarily the 
field variables within a space-time region. The 
ways differ from each other in the choice of quan- 
tities kept fixed on the 3-surface bounding this 
region. 
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It was wrong in electromagnetism to think of 
the f uy 48 fixed on the surface. It appears 
equally wrong in geometrodynamics to think of 
the g#” as fixed on the surface. 

In electromagnetism we know the potentials 
which should be fixed on the surface. In geo- 
metrodynamics we do not yet know these poten- 
tials. 

Therefore, we regard the problem ahead as 
not to change the variational principle Eq. (16), 
which is so plausible on physical grounds, but 
rather to find the right superpotential in terms 
of which to express the g#. 

In discussions (1) of the quantization of general 
relativity and (2) of the “true variables” of gen- 
eral relativity, the hope has often been expressed 
to find these superpotentials. We now have a 
third and perhaps even stronger motive to con- 
struct these quantities: to validate—if possible — 
a variational principle which, for the first time, 
summarizes in purely geometrical terms the 
whole content of classical source-free electro- 


magnetism and general relativity. 

I am deeply grateful to Professor D. C. Spen- 
cer, and especially to Professor J. A. Wheeler, 
for many long and helpful discussions on this 
problem. 
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POSSIBLE METHODS FOR THE DETERMINATION OF STRANGE- PARTICLE PARITIES* 


A. Sirlin and R. SpitzerT 
Physics Department, Columbia University, New York, New York 
(Received June 16, 1959) 


The question of the intrinsic parities of the 
strange particles has been extensively discussed 
in the past. The aim of this note is to point out 
two possible experiments that may shed light on 
this matter. 

We first consider the capture reaction 


K~ 4+d-Y4N. (1) 


It has been recently suggested that the capture 
of K” in liquid hydrogen and deuterium occurs 
predominantly from bound s states.’ In this note 
we shall consider that statement as a working 
hypothesis and shall derive a theorem that may be 
used to obtain information regarding the parity 
of the K meson relative to the YN system. With 
the assumption that K capture occurs predomi- 
nantly from s states (bound or from the continuum) 
relative to the center of mass of the deuteron, the 
most general form of the transformation matrix 
in spin space is 


M=[a+0(6,-K)(G,-K)]T for scalarK, (2a) 
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and 
Mz=([co,-K+do,-K]T for pseudoscalar K, (2b) 


where K is a unit vector along the direction of the 
relative momentum of the final system and T is 
the triplet projection operator. The operator T 
simply describes the spin correlation of the ini- 
tial state. We also note that the other possible 
invariants 6,-0,7 and 6,+(Kx6,)T can be reduced 
to the general forms given in Eqs. (2a) and (2b), 
respectively. 

We now evaluate the correlation between the 
components of the spins of the Y and N along a 
direction n, perpendicular to E: 


{0,n%n) = Tr[MM t 0,2 %n |/Tr[MM t ]. (3) 


Inserting expressions (2a) and (2b) into (3), we 
obtain 


(o,n%n) = |a - bI?/[la+ 61? +2(1al? + 101?) ]20, (4a) 
and 
(0,n%n) = - le - dl?/[le+d|?+2(Ic1*+ Id|)] <0, (4b) 
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for scalar and pseudoscalar K, respectively. 

Thus, the K parity can be determined, in prin- 
ciple, by the measurement of the sign of this 
correlation. 

The determination of (0,,0,,) in reaction (1) 
presents two main practical difficulties: in the 
first place, reaction (1) occurs only in about 1% 
of all K capture reactions; secondly, if the 
final system is An, this measurement will in- 
yolve the determination of the spin direction of 
n, which may present special problems. The 
spin direction of the A may be determined from 
the angular distribution of the decay products. 

If the final system is = p, the determination of 
the spin direction of the 2 may be difficult be- 
cause of the apparently small value of the asym- 


| metry parameter in decay. 





It is perhaps interesting to observe that if, in 
the more frequent reactions 


K +d-Y+N+z, (5) 


it turns out to be possible to isolate those events 
in which the Y and N are produced with sufficiently 
small relative momentum so that they are pre- 
dominantly in a relative s state, then a similar 
theorem would hold for their spin correlation. In 
that case, in Eqs. (2a) and (2b) the vector k 

would represent the momentum of the pion rela- 
tive to the YN system and the results of Eqs. (4a) 
and (4b) would be interchanged because of the 
pseudoscalar nature of the pion field. Reactions 
of type (5) have the advantage that they are more 
frequent than (1) and they would permit the in- 
vestigation of the spin correlations of the Ap and 


the =*n systems, which may avoid some of the 
difficulties mentioned above. 
Finally we mention the reaction 


=* + Het, He? + 2°, 


(6) 


as a possible means of determining the relative 
A-= parity. (We assume spin 3 for ,He*.) When 
the 2 has a polarization P,, the differential cross 
section for reaction (6) is given by 


do/dQ = (da/dQ), (1 + P,-P,), 


where P, is the polarization of the aHe® in the 
case P,=0 and (do/d®), is the corresponding dif- 
ferential cross section. The +(-) sign in Eq. (7) 
applies in the case in which the relative intrinsic 
parity of the initial and final systems is even 
(odd). Thus, the relative intrinsic parity of the 
final and initial systems in (6) can be determined, 
in principle, by a measurement of P, and the 
left-right asymmetry for P,#0. The derivation 
of Eq. (7) follows the same lines as the theorems 
recently given by Bilenky* and Bohr.* 

The authors wish to thank J. Franklin, N. Kroll, 
B. Margolis, and J. Steinberger for useful con- 
versations. 
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ELECTROMAGNETIC TRANSITIONS BETWEEN » MESON AND ELECTRON* 


S. WeinbergT 
Columbia University, New York, New York 


and 


G. Feinberg! 
Brookhaven National Laboratory, Upton, New York 
(Received June 15, 1959) 


The existence of the ordinary p decay, 
t~@+v+v, seems to prove that the muon and 
electron do not differ in any quantum numbers.’ 
It follows that weak electromagnetic transitions 
between muons and electrons could occur, if 
there is a mechanism to produce them. For ex- 
ample, one such mechanism would exist if the 


p. decay was not caused by a direct fevy Fermi 
interaction but instead involved a virtual charged 
boson. This particular possibility seems ruled 
out, since the predicted? rate for p ~e +y would 
be considerably greater than the upper limit set 
by recent experiments.*»* The purpose of this 
note is to discuss phenomenologically (without 
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attachment to any specific mechanism) other 
kinds of electromagnetic transitions between 
muon and electron that may be possible even if 
-e+y is somehow suppressed. 

Since muon and electron both have spin $, the 
general y.-e-y interaction may be described by 
four form factors, which correspond roughly to 
distributions of E0, MO, E1, and M1 transition 
moments. They are defined by*® 


m 
(e iJ, 0) lu) =- waa roll +a) 


q,m q" 
+f wots? + ngs) +f £1'7)"5°_xn ms 
L 


q! 
ed Oe (1) 


Here J, is the electromagnetic current, and 
9) =Pey-P pr is the momentum transfer. The 
form factors fro, f yo fe1> fyi are dimen- 
sionless functions of 
A 
P=9,9 =-F+(Q"h. 

Equation (1) follows from Lorentz invariance, 
current conservation, and the Dirac equations 
for ug and u,,, with the electron mass neglected. 

In order to avoid an infinity in (1) for ¢@-0, we 
must have 


f p90) =9, (0) =0. (2) 


to 
Furthermore, if €, is the polarization four -vector 
of a real photon with momentum -q) (¢ =0), we 
must have g,€*=0. We therefore have the famil- 
iar result that the E0 and MO terms cannot con- 
tribute to the rate of emission of a real photon, 
which is given by 


w(u -e +y) =m a8", 
B= BLU py (ODI? + IF 541 (0)1?]. (3) 


Thus the observed absence of »p ~e +7 decay tells 
us nothing about the possible existence of an E0 
or MO form factor, or about the possibility that 
fei(Y) or f y4(¢) #0, although fp 1 (0) =f y41(0) =0. 
Actually, because of the presumably “chiral” 
character of lepton interactions, it is not diffi- 
cult to find examples of theories in which we ex- 
pect that fr and f MO would be much greater 
than fry A fy1- Suppose, for instance, that in 


addition to the usual yp -decay Fermi interaction, 
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the muon and electron also have a Fermi inter- 
action of the form 


H=gy, % (1+7,)¥ udp? at +750 (4) 


with some charged particle f.° This interaction 








will generate, through “vacuum polarization” 
type diagrams, monopole form factors given by 


Fgo'P) =F yo) fl 2- vee m7) [| (5) 


where m, is the mass of f, Disa logarithmically 
pomaatier | constant, given in terms of a cutoff ) | 


by 
24)? m 
p=n(—L ; +5 i -1, (6) 
f f 


and V is’ the vacuum polarization function 


s 1 iad 1" (1+1/x)”-1 
-V(x)=5- (i a=)(1*3) of 


(1) | 








- ty as x ~-0. 
For example, if we take 7 =m)", and A=my 
then fro = -f yo =f: 
f =-g(m y/12n*) x4. 4. 


On the other hand, Ly this model fri and tun 

are of order ge’, since a virtual shite mol be 
emitted and absorbed to produce any dipole transi- 
tion moment.® This is a direct consequence of 
the “chirality conservation” of the y-e covariant 
in (4), and would still hold if the factor i 





v7 x +Y 3 


were replaced by a more complicated function of 
the f-field operators. It should be mentioned 
that, realistically, f could only be a muon ifthe ; 
coupling constant in (4) is comparable with the 
Fermi constant. This is because we must have 
ms >m p/2 in order to avoid the unobserved decay 
process p-e+f+f, and furthermore f can only 
interact weakly with nuclei in order to avoid the 
unobserved nuclear absorption process p +N-e +N’ 
in amounts comparable to ordinary p capture. 

In order to test the possible existence of the £0 
and MO terms in (1) or of the E1 and M1 terms 
for ¢ #0, it is necessary to consider cases whert 
the » ~e transition is accompanied by an inter- 
action with a virtual photon.* For example, 4 
u” bound in a 1S state around a nucleus N may 
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absorb a virtual photon from the Coulomb field 

of the nucleus, forming an electron and a re- 
coiling nucleus N’. The kinematics of this will 

be almost identical to those in ordinary p~ ab- 
sorption, the electron and N’ being emitted with 
opposite momenta of magnitude m c= 106 Mev/c, 
providing that the excitation of the nucleus is 


small. 
The rate of such a transition is given by 


- 4,5¢ 2 2 
=16m 2 op E ZI Fy | (8) 


®N-N’ 'N'? 


where Fyy»y 18 defined by 


(N' 15 °(0) IN) = (22) "*ZeF y, .(@"), (9) 


9 = all F pg (tm 7) +f yyy m7)!” 


+ IF go 2) +h py m0"). (10) 
[We are using a plane wave for the electron, a 
constant wave function for the muon, and take 
the density of states for the y~ at the nucleus 
to be Zegt*a*m,,°/nZ. Only virtual “timelike” 
photons are assumed to contribute. Terms of 
order m,,/my or (my, - my)/m, are neglected 
throughout. | 

For experimental purposes, it is best to com- 
pare (8) with the rate for ordinary y~ absorption®: 


2 
Wapg 7 (1/20°)Z ga E, mi (11) 


2_1 2 2 2 2 4 
5 =3(IC,,| +1C,,| +31C,1 +31C,,| )m . (12) 


We have then 


= 39722 2p 2/p 2 
yn “abs 3217 Z\F uy! =) /%, 


~ 2, 2 2 
=0.017Z 1 FV, .)17E, /k . (18) 
The matrix element (9) is just the Born approxi- 


mation nuclear matrix element for electron scat- 
tering (e +N ~e +N’) at momentum transfer 


(7)? =2E_sin(@/2) “me (14) 


The cross section for this process in Born approx- 
imation is 


(15) 





2 2 2 2 
(=) 2 a*cos (0/2)|F | 
* 4 
dQ lab AE sin (6/2) 
If the Born approximation gave a correct descrip- 


tion of electron scattering, one could use the ex- 


perimental data directly to compute the variation 
of Wy; _ y7/Waps from nucleus to nucleus. Since 
the Born approximation is not accurate over the 
whole periodic table, it seems more reasonable 
to use the information obtained about nuclear 
charge distributions in the electron scattering 
experiments to compute Fyn: For the elastic 
processes, where N=N’, we estimate Fyy 
using a “Fermi shape” for the nuclear charge 
density,” with radius 1.08 x1075A™“5 cm and sur- 
face thickness t=2.5x107'* cm. Then setting 
(q*)”* =5.36 x10" cm™, we obtain 


pe 3 (1.10sin6 -0.866cosé) 
NN @ 1+0.916/67 ? 


© =0.579A™. (16) 


The factor Z IFy |? in (13) then reaches a maxi- 
mum of about 6 for nuclei near Cu, falls off 
slowly for heavier elements to about 2.5 for Pb, 
where our approximations are not valid anyway, 
and falls off rapidly for lighter elements, where 
the shape probably differs from the Fermi shape 
anyway. 

In addition to the coherent process where N=N’, 
we must also consider the possibility of inelastic 
processes, where N’ is some excited state of N. 
Here we can make use of a theorem known in the 
study of muon and electron scattering."* We make 
the following assumptions: (1) Born approxima- 
tion, (2) closure approximation, (3) neglect of 
structure and mesonic effects, and (4) neglect 
of nucleon-nucleon correlations in N, except for 
the Pauli principle. Then 





2 
Dare ony! Fageny! <1/Z, (17) 
so that 
2 < 2 
Liye! Fy! Z 1+ZIF. | ‘ (18) 


This factor then reaches a maximum of about 7 
for nuclei near Cu. Since the coherent process 
N-N occurs for such nuclei at least six times 
more frequently than all inelastic processes 
(providing the above assumptions are justified), 
the electron emitted will have a very sharp 
spectrum, peaked at 105 Mev. 

It should be noted that the ordinary y™ capture, 
which involves the change of a proton into a 
neutron, cannot lead to the same nucleus, and 
thus gets no contribution from the coherent 
process. 

The best experimental evidence on possible 
electrons from y capture was obtained for cop- 
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per.” This experiment showed that 


-4 
°n - nila) / Waps(C) <5x10", 
2 
so that, taking 7 as the value of Ue! Fy! 2 
we have 


2 2 “3 
§5 [ky <4x10™. (19) 
This indicates the amount by which the up -e -y 
form factors must be reduced compared to the 
y decay coupling on the basis of present experi- 
ments. For comparison, the intermediate vector 
boson theory gives an answer of the order of 
£,”/,* ~(1/64n*)’ for this ratio. In view of the 
mysteries surrounding the relation between muons 
and electrons, a more intensive study of whether 
electrons come from yp capture would be desir- 
able. 

The experiment of Steinberger and Wolfe was 
performed to test the possible existence of yepp 
and penn Fermi interactions. The ratio of 
+N -~-e+N’ to ordinary » absorption would, for 
vector coupling of equal strength to p and n, be 
proportional to >) yA? Fy \?/Z. According to 
our calculations, this factor is about equal to 
30 for Cu. Since this enhancement factor was 
taken to be about 4 in reference 12, the sensitiv- 
ity of that experiment as a test of equal pepp 
and penn vector interactions is about 7 times 
better than quoted. 

Another experimental test of » ~e transitions 
involving virtual photons would be the internal 
conversion process p ~e+(e*+e-). The rate for 
this process depends in a complicated way on 
the form of the functions fro Suo Sev fui 
since it does not take place at fixed momentum 
transfer. (We cannot use the approximation of 
Kroll and Wada"* here, as was done in reference 
2, since we are supposing that the decay p~e ++ 
does not occur.) An order-of-magnitude estimate 
gives 


w(u~-3e) /w(u—e + +7) ~ (e*/A)(E,°/C y- 


Experimentally,"* this ratio is less than 10°, so 
that £,*/&,7<10-*. Therefore the two types of 
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experiment are about equally sensitive as tests 
of up -e -»y interactions, but the former has a 
simpler dependence on the form of the inter- 
action. 

The authors would like to thank Dr. D. Berley, 
Dr. J. Lee, Dr. R. Serber, and Dr. J. Weneser 
for helpful discussions. 
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In this section are printed the abstracts of Articles that 
have been forwarded to The American Institute of Physics 
for publication in THE PHYSICAL REVIEW. In quoting 
information obtained from this section before the appear- 
ance of the corresponding Article, reference should be 
made to “Physical Review (to be published)” rather than 
to this Journal. 


ROLE OF FUSION CHAIN REACTIONS IN THE 
NONSTATIONARY EVOLUTION OF STARS — 
SUPERNOVA STARS. M. Gryzinski, Institute 

of Nuclear Research, Polish Academy of Science, 
Warsaw, Poland (Received November 13, 1958; 
revised manuscript received May 20, 1959). 


The hypothesis is suggested that the cause of 
the explosion of supernova stars is the develop- 
ment of a fusion chain reaction in stellar material 
containing a large amount of He® at high densities, 
~2x10* g/cc. It is suggested that the conditions 
suitable for a strong concentration of He® exist 
only for relatively small hydrogen stars. The 
strength of the explosion is proportional to the 
amount of He* amassed which strongly depends 
on the mass of the star. The proto-stars with 
masses smaller than the mass of the sun lead to 
supernovae of type I and the proto-stars with 
masses equal or greater than the mass of the sun 
lead to supernovae of type II. The remnant of a 
supernova of type I is the gaseous nebula of 
luminosity due to the decay of Be” created during 
the explosion in the reaction He*(He’, 7)Be’. 


DISSIPATION OF CURRENTS IN IONIZED ME- 
DIA. O. Buneman,* Stanford University, Stan- 
ford, California (Received February 6, 1959). 


The destruction of electron drifts by instabili- 
ties is analyzed. The fastest stable drift is cal- 
culated (drift energy 0.97) and the energy of a 
faster drift is found to be dissipated into instabil- 
ities within, typically, 30 plasma periods. The 
growth of a local disturbance in this process is 
shown to take place without effective propagation. 
The “turbulent” flow pattern created, eventually, 
under nonlinear conditions is calculated numer- 
ically, demonstrating the tendency towards ran- 
domization of the initial drift energy. The effect 


stops “runaway” in about 100 plasma periods 
after which there is “heating” by “collective 
collisions” instead. 


*Now at Peterhouse, Cambridge, England. 


ELECTRON AND ION RUN-AWAY IN A FULLY 
IONIZED GAS. H. Dreicer, University of Cali- 
fornia, Los Alamos Scientific Laboratory, Los 
Alamos, New Mexico (Received February 13, 
1959). 


Hydrodynamic equations are used to describe 
the flow of the electrons and ions of a fully 
ionized gas under the action of an electric field, 
E, of arbitrary magnitude. The dynamical fric- 
tion force exerted by the electrons and ions upon 
each other through the agency of two-body Coulomb 
encounters is evaluated. In this connection the 
electrons and ions have been assigned Maxwellian 
velocity distributions which are displaced from 
each other by their relative drift velocity. This 
treatment yields a dynamical friction force which 
maximizes when the relative drift velocity is 
equal to the sum of the most probable random 
electron and ion speeds. For relative drift veloc- 
ities in excess of this value the friction force 
decreases rapidly. As a consequence, it is found 
that a fully ionized gas cannot exhibit the steady- 
state behavior characterized by time-independent 
drift velocities which has previously been ac- 
credited to it by other authors. Instead, it is 
shown that the electron and ion currents flowing 
parallel to the existing magnetic fields increase 
steadily in time (i.e., run-away) as long as a 
component of the electric field along the magnetic 
field persists. Drift velocities which greatly 
exceed the random speeds of the plasma particles 
can be created in this manner. 

The theory yields a critical electric field param- 
eter, E,, which is proportional to the plasma 
density and inversely proportional to its tempera- 
ture. It is a measure of the electric field which 
is required if the drift velocities are to increase 
and exceed the most probable random speeds in 
the gas in one mean free collision time. For 
electric fields in excess of E.,, run-away pro- 
ceeds even faster. In smaller fields, run-away 
occurs when Joule heating has depressed E.. 
sufficiently. Several interpretations of E c are 
given in terms of the collisional phenomenon 
involved. 

Within the framework of the hydrodynamic 
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equations it is shown that the well-known (tem- 
perature)” electrical conductivity law can be 
recovered, provided EXE, and the electron tem- 
perature is held constant. 

Numerical solutions giving electron tempera- 
ture and drift velocity as a function of time are 
presented for a range of the ratio E/E... The 
assumption of the displaced Maxwellian distri- 
bution is justified on the basis of a comparison 
between the rate of Joule heating and the rate of 
equipartition of random speeds. Moreover, it 
is found that the use of an anisotropic velocity 
distribution does not affect the run-away phenom- 
enon in any important way. 

The possibility of run-away induced across 
magnetic fields by steep pressure gradients and 
its relation to diffusion across magnetic fields 
is examined and discussed in detail. 


ELECTRON DAMAGE THRESHOLDS IN InSb. 

F. H. Eisen and P. W. Bickel, Atomics Interna- 
tional, Division of North American Aviation, 
Incorporated, Canoga Park, California (Received 
February 18, 1959). 


Measurements of carrier removal rate and 
isochronal recovery in electron-irradiated InSb 
indicate that displacements are produced at elec- 
tron energies as low as 240 kev. Two recovery 
stages have been found and the activation ener- 
gies for recovery determined. The conductivity 
recovery in the low-temperature stage was 
found to be first-order. 


RECOMBINATION PROPERTIES OF NICKEL IN 
GERMANIUM. G. K. Wertheim, Bell Telephone 
Laboratories, Murray Hill, New Jersey (Re- 
ceived February 6, 1959). 


A study of lifetime in germanium containing 
nickel impurities has shown that the recombina- 
tion process must be interpreted in terms of a 
multilevel model, with three charge states as- 
signed to nickel. At room temperature the elec- 
tron-capture cross section of neutral nickel is 
0.96 x10~** cm?, and that of singly negative nickel 
is 5.9x107** cm?; the capture constants in both 
cases are independent of temperature. The cross 
section of doubly negative nickel for holes is 
approximately 1x10~"* cm?; no information on 
the temperature dependence of this cross section 
has been obtained. Lifetime data are in good 
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agreement with the known energy levels and goliq 
solubility. The results are compared with thoge 
in the literature. 





PARAMAGNETIC RELAXATION IN NICKEL 
FLUOSILICATE. K. D. Bowers and W. B. Mims, 
Bell Telephone Laboratories, Murray Hill, New 
Jersey (Received February 17, 1959). 





Physical models of the paramagnetic relaxa- 
tion process are briefly reviewed, and a descrip. 
tion is given of the pulsed-microwave method of 
investigation. This method has been applied to 
the study at low temperatures of zinc/nickel 
fluosilicate at various magnetic dilutions, and 
an additional relaxation process has been ob- 
served which is associated with the sharing of 
energy among the spins. 

The results depend only on the temperature 
and magnetic concentration, and are independent | 
of all other experimental conditions provided that 
excessive power levels are avoided. The lack of 
size dependence indicates that the relaxation 
rate is not controlled by spatial diffusion of 
phonons. The temperature dependence which has 
been observed is too acute to fit the 1/T law of 
the simple “direct process” model, and suggests 
that the relaxation rate is controlled either by 
phonon combination probabilities or by a Raman 
type of interaction between phonons and spins. 





ELASTIC CONSTANTS OF a-BRASSES: VARIA- | 
TION WITH SOLUTE CONCENTRATION FROM | 
4.2 TO 300°K. J. A. Rayne, Westinghouse Electric, 
Corporation, Pittsburgh, Pennsylvania (Received | 
February 6, 1959). 


Elastic-constant measurements have been made | 
on a representative series of a-brasses in the 
temperature range 4.2-300°K. No appreciable 
change with temperature occurs in the dependence 
of C and C’ on solute concentration. Possible 
explanations for this result are considered. The 
€ values for the alloys, as computed for the 
elastic-constant data at 4.2°K, differ significantly 
from the calorimetric values. It is believed that 
this discrepancy results from the difference be- 
tween the phonon wavelengths used in the ultra- 
sonic measurements and those contributing to 
the lattice heat capacity of the alloys at low 
temperatures. 
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ANTIFERROMAGNETIC STRUCTURES OF MnS,, 
MnSe,, and MnTe,. J. M. Hastings, N. Elliott, 

and L. M. Corliss, Chemistry Department, Brook- 
haven National Laboratory, Upton, New York 
(Received February 13, 1959). 


The antiferromagnetic structures of MnS,, 
MnSe,, and MnTe, have been obtained by neutron 
diffraction. The disulfide exhibits ordering of 
the “third” kind, the ditelluride, ordering of the 
“tirst” kind, and the diselenide, an arrangement 
which is intermediate between the two. The 
structures are discussed from the viewpoint of 
indirect exchange. Magnetic susceptibilities, ob- 
tained by the Gouy method, are reported for all 
three compounds. 


ELECTRON PARAMAGNETIC RESONANCE IN 
$rTiO,. W. L Dobrov, R. F. Vieth, and M. E. 
Browne, Lockheed Research Laboratories, Palo 
Alto, California (Received February 13, 1959). 


The intensities of paramagnetic resonance 
lines in SrTiO, have been correlated to the iron 
concentrations. Evidence for attributing the 
origin of the resonance in titanates to iron im- 
purities (Fe**) rather than to the domain struc- 
ture was obtained. In SrTiO,, resonance ex- 
periments over the temperature range 300°K to 
1.9°K show a slow continuous increase of the 
tetragonal component of the crystalline field 
down to the ferroelectric region. From intensity 
considerations the absolute signs of the crystal- 
line splitting parameters have been determined. 


ENERGY DEPENDENCE OF RECOVERY IN 
IRRADIATED COPPER. J. W. Corbett and R. M. 
Walker, General Electric Research Laboratory, 
Schenectady, New York (Received February 9, 


1959). 


Stage I recovery (10°K-65°K) in electron-irradi- 
ated copper consists of five substages of recovery, 
designated as I (4 —Ig in order of increasing tem- 
perature. We report here the shift in the popula- 
tions of I [4 Ig as a function of the energy of the 
bombarding electrons. The sum of the first 
three substages (I [4> Ip, and Ic) increases with 
decreasing bombarding energy at the expense of 
the higher temperature substages. The ratio 
D/A, also increases with decreasing energy. 

The data are interpreted as corroboration of a 


previously proposed model for stage I recovery. 
The measured shifts are quite small. This is 
discussed in relation to the recoil energy dis- 
tributions and it is concluded that the most likely 
interpretation is that the average separation be- 
tween a primary recoil and its vacancy does not — 
vary rapidly with increasing recoil energy near 
threshold. Several possible explanations for 

this behavior are outlined. 


X-RAY MEASUREMENT OF THE ATOMIC SCAT- 
TERING FACTOR OF IRON. Boris W. Batterman,* 
Bell Telephone Laboratories, Murray Hill, New 
Jersey (Received February 11, 1959). 


The x-ray integrated intensities of Bragg re- 
flections of iron, copper, and NaCl powders were 
measured with monochromatic Fe K,, radiation. 
By investigating the effect of cold work on inten- 
sities of the metal powders and comparing the 
measured intensities of different orders with 
theoretical values, it was shown that extinction, 
surface roughness, and preferred orientation 
effects were negligible. The four measured iron 
reflections gave atomic scattering factors rela- 
tive to either copper or NaCl which agreed, with- 
in experimental error, with values calculated from 
wave functions of the isolated iron atom. It is 
estimated that the experimental form factor for 
the lowest order reflection of copper and iron 
agrees with free atom values to an error cor- 
responding to the contribution of + one free-atom- 
3d electron. This is to be contrasted with results 
of Weiss and DeMarco, who concluded from single- 
crystal measurements that metallic iron has 
(2.3+0.3) 3d electrons, compared with 6 for the 
free atom. 


“Present address: IBM Research Laboratory, 
Poughkeepsie, New York. 


K ABSORPTION EDGE OF SELENIUM. Richard 
D. Deslattes* and Hillery S. deBen, Florida State 
University, Tallahassee, Florida (Received 
February 9, 1959). 


Use of thin unsupported selenium absorbers 
has permitted observation of significant struc- 
ture in the K absorption spectrum. Observations 
of the Kf, emission and the fundamental optical 
absorption of similar samples have been made 
and used to assist in the interpretation of the 
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K-edge structure. It is suggested that the prom- 
inent absorption maximum at the experimentally 
determined “edge” arises from transitions to an 


x-ray excitation state or states in the region of 


the forbidden energy gap. The following general 
absorption is attributed to transitions to the con- 


duction band. The gap energy required in this 


interpretation is approximately 4 electron volts, 


indicating that the lower part of the conduction 
band contains little p-type symmetry. 


* 
Present address: Cornell University, Ithaca, New 
York. 


MAGNETIC SUSCEPTIBILITY OF CLOSE-PACKED 


HEXAGONAL GOLD-INDIUM ALLOYS. T. B. 


Massalski, * Department of Physical Metallurgy, 


The University, Birmingham, England, and 

L. Meyer and D. Weiner, Institute for the Study 
of Metals, The University of Chicago, Chicago, 
Illinois (Received February 16, 1959). 


The magnetic susceptibility of close-packed 


hexagonal gold-indium alloys has been measured 


at 300°K, 77°K, 4.2°K, and 1.2°K. The suscep- 
tibility data show a striking similarity to the 
c-axis spacing exhibiting a maximum at about 
19 atom % In, consistent with a model assuming 


that the Fermi surfaces touch the Brillouin zone 


face (0002) at this concentration. 


a 
Now at Mellon Institute, Pittsburgh, Pennsylvania. 


APPROACH TO BOLTZMANN EQUILIBRIUM OF 


THE NUCLEAR SPINS IN A QUADRUPOLAR 
SOLID. Walter I. Goldburg, Carnegie Institute 


of Technology, Pittsburgh, Pennsylvania (Received 


August 21, 1958; revised manuscript received 
April 1, 1959). 


A single crystal of NaNO, is oriented in an 


external magnetic field at an angle (6,) for which 


the quadrupole splitting is large. With radio- 
frequency radiation, a pair of Zeeman levels of 


the sodium nuclear spin system (J = 3) is saturated. 


Then the crystal is quickly rotated to an angle 
where the inequality of the level spacings is a 
minimum (approximately 4 the line width) and 
where spin mixing can take place. The crystal 
is then quickly returned to @, where, by nuclear 
resonance pulse techniques, the populations of 
the four sodium levels can be determined. This 
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id 
sequence of operations is completed in a time 
short compared to the spin lattice relaxation time 
When the m = 3 levels are initially saturated, 
the spin system rapidly comes to the expected 
spin temperature, 7, =4f the lattice temperature. 
When the outer satellites are saturated, however, 









the resulting equilibrium populations of the levels 


are not those of a Boltzmann distribution. Algo 
measured was the rate of approach of the spin 
system to equilibrium for various larger split- 
tings. For splittings of the order of the line Width, 
an unexpected plateau appears in the spin mixing 
rate. The usual rate equations do not seem ade- 
quate to describe the observed dynamic behavior 
of this spin system. 


HARTREE-FOCK THEORY WITH NONORTHOG- 
ONAL BASIS FUNCTIONS. R. McWeeny, Depart- 
ments of Mathematics, Physics and Chemistry, 
University College of North Staffordshire, Staf- 
fordshire, England (Received January 19, 1959), 


Solutions of the Hartree-Fock self-consistent- 
field equations can be approximated by linear 
combination of a set of basis functions: the equa- 
tions then assume a matrix form, as shown by 
Roothaan. It is possible, however, to obtain the 
required variational solution of the many-elec- 
tron problem much more directly by an iterative 
construction of the density matrix. This method, 
first developed by the author for expansion in 
terms of orthogonal basis functions, is extended 
to the case of a nonorthogonal set. 


GROUND STATE OF THE HELIUM ATOM. II. 
Toichiro Kinoshita,” Bell Telephone Laboratories, 
Murray Hill, New Jersey, and Laboratory of 
Nuclear Studies, Cornell University, Ithaca, 
New York (Received February 19, 1959). 


A further attempt is made to improve the the- 
oretical prediction of the energy of the ground 
state of atomic helium. The nonrelativistic part 
is treated by the variational method of Stevenson 
and Crawford which is useful for improving the 
lower bound for the ground state energy. Linear 
combinations of up to 80 terms of generalized 
Hylleraas type are employed in the numerical 
computation. The best trial function gives 
-2.9037237 as an upper bound and -2.9037467 as 
a lower bound for the ground state energy. It is 
estimated from the calculated results that the 
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exact nonrelativistic energy of He ground state 
will be found in the neighborhood of -2.9037247. 
Rigorous formulas are derived which can be used 
for calculating the upper limits to the errors in 
the expectation values of mass polarization and 
relativistic corrections. Although these formulas 
give very broad limits of error, they are useful 
in estimating the order of magnitude of actual 
errors ina semiempirical manner. With mass 
polarization and relativistic corrections as well 
as electrodynamical corrections, the theoretical 
ionization potential becomes 198 310.76 cm™ 

which is in good agreement with the latest ob- 
served value 198 310.8, +0.15 cm™. 


* 
Permanent address: Cornell University, Ithaca, 
New York. 


| TOK ELECTRON CAPTURE RATIO IN Ge”. 
D. E. Rehfuss* and Bernd Crasemann, Depart- 
ment of Physics, University of Oregon, Eugene, 


| Oregon (Received July 25, 1958). 


By comparison of the number of 0.175- Mev 
gamma rays occurring in the decay of As” with 
the number of gallium K x-rays proceeding from 
Ge” K capture, the relative amount of L capture 
inGe™ has been determined indirectly. Solid 
sources and a scintillation spectrometer were 
employed. The Ge™ 1/K electron capture ratio 


' is found to be 0.09+0.05, in agreement with 
| theory. 


"Now at the Los Alamos Scientific Laboratory of the 
University of California, Los Alamos, New Mexico. 


NEUTRON- PROTON CAPTURE CROSS SECTION. 
Robert W. Stooksberry and Marshall F. Crouch, 
Case Institute of Technology, Cleveland, Ohio 
(Received January 28, 1959). 


The mean lifetime of thermal neutrons in water 
is measured with a large moderator (radius ~7 
diffusion lengths), for which the correction for 
escape is only about 5%. The geometry is such 
that the perturbing effect of the immersed BF, 
neutron detector can be almost rigorously cal- 
culated. The value obtained is 7=206.3+5.0 
usec, giving the neutron-proton capture cross 
section as (0.330 + 0.008) x10-*4 cm?. A possible 
systematic error from higher modes in the neu- 
tron distribution is discussed. 





INVESTIGATIONS OF THE Be*(p, )B® and Be*(, 
ay)Li® REACTIONS. Jerry B. Marion, Los 
Alamos Scientific Laboratory, Los Alamos, New 
Mexico, and University of Maryland, College 
Park, Maryland, and Jules S. Levin, Los Alamos 
Scientific Laboratory, Los Alamos, New Mexico 
(Received February 12, 1959). 


A pulsed-beam time-of-flight technique has 
been used to study neutrons and y rays from the 
proton bombardment of Be® in the energy range 
from 2.3 to 5.4 Mev. Excitation curves for the 
ground-state neutron group from the Be*(p, n)B® 
reaction have been measured at 0° and 90°. 
These curves show the well-known resonances 
at 2.56 and 4.6 Mev and, in addition, show a 
broad maximum near 3.5 Mev. The y-ray exci- 
tation curve also indicates maxima at 2.56, 3.5, 
and 4.5 Mev. The 3.5- and 4.6-Mev resonances 
are tentatively assigned to p- and d-wave states 
having (J =2+, T=1) and J/=3', T=1), respec- 
tively. The previously observed maximum near 
4.9 Mev appears to be the result of the effect of 
the 3.5-Mev resonance on the n, group immedi- 
ately above its threshold. The cross section for 
the production of continuum neutrons is appre- 
ciable for proton energies above about 3 Mev. 

A large fraction of these neutrons appears to re- 
sult from the neutron decay of the 2.43-Mev Be® 
level following the inelastic excitation of this 
level. 


NUCLEAR SPECTROSCOPIC STUDIES OF 
NEUTRON-DEFICIENT ISOTOPES IN THE RARE 
EARTH REGION. Kenneth S. Toth* and John O. 
Rasmussen, Radiation Laboratory and Depart- 
ment of Chemistry, University of California, 
Berkeley, California (Received November 3, 
1958). 


The conversion-electron spectra of four neu- 
tron-deficient rare earth isotopes are presented 
and discussed. The four isotopes are Tb’*’, 
Tb'**, Dy'®®, and Ho’. Mihelich, Harmatz, and 
Handley in a recent publication have presented 
the conversion-electron data of a large number 
of rare earth isotopes, including all but one 
(Dy'**) of the four above-mentioned activities. 
The agreement between the two works is good. 
Some additional conversion lines, not reported 
by Mihelich et al., have been identified in this 
work as belonging to the decay of Tb'**, Tb’, 
and Ho’®’. Decay schemes are proposed for 
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Dy"* and Tb"*. Tb*** and Ho’ decay to levels 
in Gd'* and Dy’®, respectively. The electron- 
capture decay data of Tb'** and Ho'™ are com- 
pared with data obtained by other workers in the 
negatron decay of Eu'™ and Tb’® to levels in 
Gd'™ and Dy'®, respectively, and transitions 
common to both modes of decay are noted. 


"Present address: University Institute for Theoreti- 
cal Physics, Copenhagen, Denmark. 


NEW TERBIUM ISOTOPE, Tb***. Kenneth S. 
Toth,* Kenneth T. Faler, and John O. Rasmussen, 
Radiation Laboratory and Department of Chem- 
istry, University of California, Berkeley, Cali- 
fornia (Received November 6, 1958). 


A careful examination of photon and electron 
spectra from cyclotron-produced mixtures of 
light terbium isotopes has led to the indentification 
of a new isotope, Tb'*, with an 18.5-hr half-life. 
The 344.1 -kev transition from first excited state 
to ground in daughter Gd'*? has been identified 
and shown to be identical in energy to a tran- 
sition associated with the beta decay of 9.3-hr 
Eu'**™ to the same daughter nucleus. Gamma 
transitions of 180 and 265 kev are also assigned 
to the decay of Tb'**. 


* 
Present address: University Institute for Theoreti- 
cal Physics, Copenhagen, Denmark. 


BETA DECAY OF Bk**’ AND Bk***. Susanne E. 
Vandenbosch, Herbert Diamond, Ruth K. Sjoblom, 
and Paul R. Fields, Argonne National Laboratory, 
Lemont, Illinois (Received February 4, 1959). 


The beta decay of Bk*® has been studied with a 
beta spectrometer. Bk**° has two beta groups: 
725415 kev (89+1% abundant) and 1760+ 50 kev 
(11+1% abundant). Conversion electron lines 
corresponding to 42.2, 98.2, 890, 930, 990, and 
1032-kev transitions were observed and their 
intensities measured. The information obtained 
from crystal spectrometer singles and coinci- 
dence measurements was combined with beta 
spectrometer data to construct a decay scheme 
for Bk*** involving gamma vibrational levels in 
Cf**° analogous to those found in Pu***. The beta 
half-life of Bk*®° is 193.3+0.3 minutes. The 
electron-capture partial half-life of Bk’™ is esti- 
mated to be greater than 5 years. The beta 
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ee, 


spectrum of Bk**® shows an allowed or first for. | 


bidden transition with an end-point energy of 
125+2 kev. 


METASTABILITY OF PLATINUM-199. Morris 4 
Wahlgren and W. Wayne Meinke, Department of 
Chemistry, University of Michigan, Ann Arbor, 
Michigan (Received February 12, 1959). 


A metastable state of platinum-199 has been 
produced by thermal neutron irradiation of norm) 
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life 


and enriched platinum samples. The isomer de. 


cays with a half-life of 14.1+0.3 seconds by the 
emission of y rays of 32+ 2 and 393 +2 kev ener. 
gy. The thermal neutron activation cross section 
of Pt’ for the formation of the isomer is 0.028 
+0.003 barn. Tentative level assignments are 


made, consistent with systematics and shell —_ 


ANGULAR DISTRIBUTIONS OF NEUTRONS IN- 
ELASTICALLY SCATTERED FROM SODIUM. 
E. N. Shipley,* G. E. Owen, and L. Madansky, 
Physics Department, The Johns Hopkins Univer- 
sity, Baltimore, Maryland (Received February 
6, 1959). 


The angular distribution and total cross sec- 
tion for inelastic neutron scattering proceeding 
to the first excited state in sodium have been 
measured by observing the neutron-gamma coin- 
cidences. The measured total cross sections 
for the various neutron energies used are 
o (3.49 Mev) =370+120 mb, o7(3.75 Mev) =270 
+90 mb, and o7(4.00 Mev) =280+90 mb. The 
differential cross sections were found to be 
higher in the backward than in the forward direc- 
tions. 


* 
Now at Physics Department, Northwestern Univer- 
sity, Evanston, Illinois. 


DECAY OF SHORT-LIVED BARIUM AND LAN- 
THANUM FISSION PRODUCTS. R. P. Schuman, 
E. H. Turk, and R. L. Heath, Phillips Petroleum 
Company, Atomic Energy Division, Idaho Falls, 
Idaho (Received February 13, 1959). 


The half-lives and gamma spectra of the short: 
lived barium and lanthanum fission products 
have been measured. Ba’ has a half-life of 18 
+1 min and emits x-rays and gamma rays of 
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120 kev, 190 kev, 290 kev, ~350 kev, 460 kev, 
t (40 kev, and 740 kev. Ba’ has a half-life of 
{1+1 min and emits x-rays and gamma rays of 
90 kev, ~260 kev, 890 kev, 970 kev, 1080 kev, 
1200 kev, 1860 kev, and 1680 kev. La’ has a 
half-life of $.85+0.10 hr and emits a low-inten- 
is A sity gamma ray at 1370 kev. La’* has a half- 

a | life of 85+ 6 min and emits gamma rays of 640, 
oT, § 990, 1080, 1080, 1370, 1540, 1750, 1920, 2080, 
2400, 2570, 3000, 3300, and 3650 kev. 


- 


~ 


ENERGY RESONANCES. P. A. Egelstaff, Atomic 
Energy of Canada Limited, Chalk River, Ontario, 
tion Canada (Received February 9, 1959). 


de- | STATISTICS OF SLOW-NEUTRON “NEGATIVE” 


2 Neutron cross sections measured with very 
hy slow neutrons often reveal the influence of levels 
" below the neutron binding energy. The published 

data on 20 isotopes are analyzed in such a way 
that the properties of the levels below and above 

N- | the binding energy may be compared. It is found, 

in particular, that the neutron width and level 

’ | spacing distributions for these (- and+) levels 

ver | are the same within the limits of the present 

YF analysis. 


— 


mg | DECAY OF K®. N. Benczer-Koller, A. Schwarzs- 
child, and C. S. Wu, Columbia University, New 
oin- | York, New York (Received January 2, 1959). 


The decay of K® has been reinvestigated in an 
' attempt to clarify the discrepancies existing in 


70 | 

| the ordering of the energy levels of Ca®. The 

: external conversion spectrum showed lines cor- 
~~ responding to the following gamma-ray energies: 


_ 0.220+0.002, 0.371 40.002, 0.388+0.002, 0.394 
+0.002, 0.591+0.003, 0.614+0.004 Mev. In 
addition a y ray of 1.005 Mev was observed in 
the scintillation spectrum. 8-y coincidence 
measurements establish that the gamma transi- 
tions to the ground state are the 0.371-Mev and 
the 0.591-Mev gamma rays and not the 0.614- 

N- |} Mev gamma ray as believed previously. These 
an, | results and the y-y coincidence measurements 
cul | indicate only one possible level scheme in Ca*: 
lis, | 0, 0.371, 0.591, 0.985, 1.373 Mev. These en- 
ergy levels as well as their proposed spin and 
ort: | Parity assignments are in good agreement with 
the results from the MIT nuclear reaction data. 
18 | Atwo-cycle baffle for the solenoid 8 spectrom- 
eter is described. This baffle reduces signifi- 


er- | 





cantly the background due to scattered gamma 
radiation from very strong external conversion 
sources. 


RANGE -ENERGY RELATIONS FOR PROTONS 
IN Be, C, Al, Cu, Pb, AND AIR. R. M. Stern- 
heimer, Brookhaven National Laboratory, Upton, 
New York (Received February 6, 1959). 


Range-energy relations for protons have been 
obtained for six substances (Be, C, Al, Cu, Pb, 
and air). The calculations for the energy loss 
dE /dx include the shell corrections at low ener- 
gies, and the density effect which becomes im- 
portant in the high-energy region. The present 
results can also be used to determine the range 
of » mesons up to ~10 Bev. Besides the calcu- 
lated values of the range, tables of the ionization 
loss dE /dx are also presented. 


INTERACTION OF YTTRIUM WITH PROTONS 

OF ENERGY BETWEEN 60 AND 240 Mev. Albert 
A. Caretto, Jr.,* and Edwin O. Wiig, Department 
of Chemistry, University of Rochester, Rochester, 
New York (Received October 24, 1957; revised 
manuscript received March 30, 1959). 


Absolute cross sections are reported for a 
number of nuclides produced by the interaction 
with yttrium of protons of 60, 100, 150, 180, 
and 240 Mev. At low energies the yields can be 
accounted for by direct interaction or knock-on 
processes. At higher energies the knock-on 
cascade model together with evaporation appears 
to explain the observed yields. 


Present address: Lawrence Radiation Laboratory, 
Livermore, California. 


y-y ANGULAR CORRELATION MEASUREMENTS 
OF THERMAL NEUTRON CAPTURE y-RAYS IN 
MAGNESIUM, SILICON, PHOSPHORUS, AND 
SULFUR. G. Manning * and G. A. Bartholomew, 
Atomic Energy of Canada Limited, Chalk River, 
Ontario, Canada (Received February 16, 1959). 


A description is given of an angular correlation 
apparatus, using two Nal(T1) scintillation spec - 
trometers, which has been used to study cas- 
cades of y-radiation following thermal neutron 
capture. The elements studied were magnesium, 


121 








VoLuME 3, NUMBER 2 


PHYSICAL REVIEW LETTERS 



















JULY 15, 1959 F vou 

i — 

silicon, phosphorus, and sulfur. The following ential cross sections are interpreted with the aij Di 
new spin assignments were made: the 4.93-Mev of predictions of the Butler direct interaction pour 
level of Si?®, 3/2; 6.38-Mev level of Si?*, 1/2; (surface) theory. It is concluded that this theory | cula 


1.15-Mev level of P™, 1; and 3.26-Mev level of 
p*, 2. In addition, the spins of the 3.41-Mev 
level of Mg”® and the 3.22-Mev level of S* have 
been confirmed to be 3/2. Other aspects of the 
decay schemes are discussed. 


* 
Present address: Norman Bridge Laboratory of 
Physics, California Institute of Technology, Pasadena, 

California. 


Li’(n, t)He* REACTION AT INTERMEDIATE EN- 
ERGIES. S. J. Bame, Jr., and R. L. Cubitt, Los 
Alamos Scientific Laboratory, University of Cali- 
fornia, Los Alamos, New Mexico (Received 
November 6, 1958). 


The Li®(n, t)He* reaction cross section has been 
measured in a neutron energy interval extending 
from 9 kev to 340 kev. The cross section descends 
from 1.77 barns at 9 kev to a minimum of 0.74 b at 
100 kev before rising to a peak resonance value 
of 2.75 barns at 258 kev. The cross section 
from 9 kev to 90 kev is not proportional to E-”?, 
but can be represented by o =3.96E °°’ with o 
in barns and the neutron energy E in kev. Angu- 
lar distributions of the tritons are presented for 
neutron energies of 150, 200, 258, 300, 350, 
and 565 kev. These distributions, measured in 
a relative manner, are placed on an absolute 
scale by normalizing to the cross-section data 
of this experiment, extended to 565 kev by use 
of previously measured cross sections. 


DIFFERENTIAL CROSS SECTIONS OF SOME 
(a,p) REACTIONS. C. E. Hunting” and N. S. 
Wall, Department of Physics and Laboratory 
for Nuclear Science, Massachusetts Institute of 
Technology, Cambridge, Massachusetts (Re- 
ceived February 16, 1959). 


Using 30.4-Mev alpha particles, proton differ - 
ential cross sections have been obtained for the 
following transitions, corresponding to discrete 
states of the residual nuclei: B’(a, pe ee 
Na*(a, p)Mg”* ond, 1.83 Mev? P™(a, b)S™* on 
and Al*"(a, p)Si*gnd, 2.24 Mev- A eten — 
spectrum corresponding to the first 12 Mev of 
excitation of the residual nucleus was obtained 
for the case involving the Al”’ target. The differ- 
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is probably a more useful tool in interpreting th, ' for: 


results obtained than is the compound nucleus _ tion 


theory. Some deviations from the predictions of | eval 
the surface direct-interaction calculations are rate 
discussed. a 
* ; N 
Now at the Oak Ridge National Laboratory, Oak ) and J 
Ridge, Tennessee. San I 


; 


ENERGY LEVELS IN Na* FROM THE Ne™(p,))x GEO 
REACTION. J. J. Singh, V. W. Davis,* andR.w.} MAG 











Krone, Department of Physics, University of Mar} 

Kansas, Lawrence, Kansas (Received January 2, 
1959). = 
| first 
The gamma-ray excitation function for the certa 
Ne”*(p,y) reaction has been investigated in the gravi 
energy range of 600 kev to 1800 kev. Eleven res- Maxv 
onances could definitely be identified as resulting| in te1 
from this reaction. The cascades leading to the | are p 
ground state of Na”* were determined and found tp| trom: 
proceed for the most part through well-known vecto 
lower lying states, although no gamma rays cor- § time 
responding to transitions through the 2.08-Mev fF netic 
and 2.39-Mev states were observed at any of the | apply 
resonances investigated. or the 
* has y 
Now at Los Alamos Scientific Laboratory, Los corre 
Alamos, New Mexico. to cal 
spino: 

DISPERSIVE EFFECTS IN ELECTRON-NUCLEIS| 
INTERACTIONS. Nicholas A. Krall” andE.E. | pypy 
Salpeter, Laboratory of Nuclear Studies, Cornell ppr, 
University, Ithaca, New York (Received borate 
February 16, 1959).  Califc 
Elastic scattering of electrons from nucleiis | The 
influenced by the possibility of virtual transitions gjgoy, 





to excited nuclear states in intermediate states. 
Such dispersive corrections to elastic electron- 
deuteron scattering are calculated in second orde} ryje p 
Born approximation for incident electron energits} igoyec 
from 200 Mev to 500 Mev for various values of i that 5 
momentum transfer. The static second order Borg ang B, 
formulas are also evaluated. Similar, but less 
accurate, calculations are carried out for He’, 
C™, and some heavier elements. These results 
are used to find small corrections to nuclear radii 
obtained from an analysis of the Stanford experi 
ments. 





















2. 


~ 


5 
W. 


7 26 


eR 





EUS! 


nell 





VotuME 3, NUMBER 2 


PHYSICAL REVIEW LETTERS 





Juty 15, 1959 








Dispersive effects on the energy-level shifts in 
pound s-states of hydrogenic atoms are also cal- 
culated, using second order perturbation theory 


} ‘or the deviation of the nuclear charge distribu- 


tion from a point charge. These small shifts are 


evaluated for deuterium and He and, less accu- 


rately, for other nuclei with Z <50. 


*Now at the John Jay Hopkins Laboratory for Pure 
and Applied Science, General Atomic Corporation, 
San Diego, California. 


GEOMETRY OF GRAVITATION AND ELECTRO- 
MAGNETISM. Louis Witten, RIAS, Baltimore, 
Maryland (Received January 21, 1959). 


An independent derivation is given of equations 
first derived by Rainich which show how, under 
certain circumstances, the combined theory of 
gravitation and electromagnetism of Einstein and 
Maxwell can be unified and described exclusively 
interms of geometry. Some algebraic relations 
are presented between the Ricci tensor, the elec- 
tromagnetic field tensor, and their principal null 
vectors. It is shown that in regions of space- 
time where the two invariants of the electromag- 
netic field both vanish, the unified theory cannot 
apply. Either such regions do not exist in nature 
or their description in terms of pure geometry 
has yet to be found. Advantage is taken of the 
correspondences between tensors and spinors 
to carry out most of the present calculations in 
spinor space. 


PHENOMENOLOGICAL ANALYSIS OF HYPERON 
DECAY. S. A. Bludman, Lawrence Radiation La- 


_ boratory, University of California, Berkeley, 


rder’ 
gies 





California (Received February 16, 1959). 


The AJ=3 rule is incorporated into a previously 
discussed universal weak Yukawa interaction by 
assuming the simplest relation between chirality 
and charge operators. Expressing the Al = 4 
tule by forming the ‘ and = isospinors into an 
isovector B and an isoscalar B,, we assume 
that B occurs (as in B decay) with gy, (+3775) 
and B, (which has no B-decay counterpart) with 





$ry,). Here g is the constant previously 


BY (3 - 
fitted to the 7 decay rate, and ¢ the ratio of 
Gamow-Teller and Fermi coupling constants. 
Depending on the sign taken between the B and 
B, interaction terms, = decays into n+7n* in 


pure S and into p+7° in pure P channels, or 
vice versa. In either case, £*+~p+7° involves 
maximal S-P interference and a°=0.98. Decay 
into J= 3 proceeds via y,($+ $7y5) or y,($+ $775), 
depending on whether > decay is pure S or pure 
P. The second case, but not the first, leads to 
a A decay rate in agreement with experiment. 

In this case a, =0.54; and in = decay, a==0.64 
and the calculated decay rate is 2.4107" sec". 


INTERFERENCE EFFECTS IN LEPTONIC DE- 
CAYS. Steven Weinberg, Columbia University, 
New York, New York (Received February 17, 
1959). 


It is proven that in a leptonic decay experiment 
in which the lepton masses and charges may be 
neglected, and in which no pseudoscalar corre- 
lations are measured, all VA interference terms 
will be antisymmetric under exchange of the two 
leptons, while the pure Vand A terms will be 
symmetric. If the experiment measures a pseu- 
doscalar correlation, these conclusions are re- 
versed. Even if the lepton masses cannot be ig- 
nored (e.g., for A°~u" ++), or low-energy 
decay), it is still true that no VA interference 
may appear when scalars are measured, and 
only VA interference may contribute when pseu- 
doscalars are measured, providing that the lepton 
spins and momenta are not directly observed. 
Thus experiments can be devised that involve no 
interference effects, or only interference effects. 
This theorem holds independently of the strangeness 
change, spin change, energy transfer, or of any 
particular assumptions about the form of the V 
and A currents. It proves most useful when it 
is difficult or tedious to calculate transition rates 


directly. Applications are discussed, including 
possible tests of the Feynman-Gell-Mann theory 
in nonunique forbidden 8 decay, of the nature of 
the leptonic A° and K° decay interaction, and of 

the charge symmetry properties of weak inter - 

actions. 





